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Abstract 



A search for slepton production via the decay of pair-produced neutralinos has 
been performed under the assumption that the sleptons have observable lifetime in 
the detector before each decays to a lepton and a gravitino. Sleptons, neutralinos 
and gravitinos are particles predicted by the theory of supersymmetry, and are 
the super symmetric partners of the Standard Model leptons, neutral bosons and 
of the graviton respectively. The search was performed in 628 pb _1 of data taken 
by the ALEPH detector at LEP centre-of-mass energies from 189 to 208 GeV. It 
was motivated by general predictions of Gauge-Mediated Supersymmetry Break- 
ing (GMSB) models in which the lightest supersymmetric particle is always the 
gravitino. No evidence of the process was found. Model-independent cross-section 
limits are quoted as a function of neutralino mass, slepton mass and slepton life- 
time in the case that the neutralino branching ratios to each slepton are equal at | 
(the so-called slepton co-NLSP scenario, where NLSP stands for 'Next-to-Lightest 
Supersymmetric Particle') and in the case that the neutralino decays exclusively to 
the stau (the stau-NLSP scenario). Excluded regions in the neutralino-stau mass 
plane are shown for four gravitino masses under model-specific assumptions. 
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Chapter 1 
Introduction 



This thesis documents a search for evidence of supersymmetry, which is a theory 
relating fermions and bosons. Supersymmetry predicts that every Standard Model 
particle should have a partner which is identical in every way apart from its spin, 
which should differ by |. Thus, under supersymmetry, the existing particle spec- 
trum must be at least doubled in size to accommodate the new supersymmetric 
particles. There are many motivations for taking supersymmetry as a serious can- 
didate for physics beyond the Standard Model, but a complete and viable model of 
supersymmetry must include a mechanism for its breaking, since if the supersym- 
metric particles had the same masses as their Standard Model counterparts they 
would have been discovered some time ago. GMSB (Gauge-Mediated Supersym- 
metry Breaking) is such a model. A key prediction that makes its phenomenology 
different from other models is that the lightest supersymmetric particle should be 
the 'gravitino' (the supersymmetric partner of the graviton - the quantum of the 
gravitational interaction), and that it is very light. 

This thesis documents a search for the production of sleptons (the supersymmet- 
ric partners of leptons) via the decay of neutralinos (the supersymmetric partners 
of the Standard Model neutral bosons) that in turn are pair-produced in electron- 
positron collisions. It is assumed that the sleptons travel an observable distance 
in the detector before decaying to their Standard Model counterparts and a grav- 
itino. The primary distinguishing feature of these events in a detector would then 
be the presence of particles whose trajectories do not originate at the primary 
e + e~ interaction point. This process can have an advantage over direct slepton 
pair-production for the prospect of discovering supersymmetry since GMSB mod- 
els can predict the neutralino production cross-section to be higher than that of 
the slepton. 

Chapter ^ will introduce the theory of supersymmetry in more detail, begin- 
ning with how it can provide answers to some of the unanswered questions of the 
Standard Model, and ending with how GMSB could manifest itself in an e + e~ col- 
lider. Chapter ^| describes the ALEPH detector. Chapter f| discusses the simulated 
data sets used to develop the search techniques, and the real data set in which 
the search was performed. Chapter || will describe the techniques used to analyse 
particle trajectories found not to originate at the primary interaction point in order 
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CHAPTER 1. INTRODUCTION 



to test whether they conform either to the hypothesis that they originated from 
certain known background processes, or from slepton decay. Chapter [| gives the 
cut-based selections on event variables that are used as the final discrimination 
between signal and background. Finally, Chapter ^ gives the search outcome (that 
no evidence for the process was found), and then proceeds to discuss methods for 
quantifying this result in terms of supersymmetric parameter space, before giving 
model-independent cross-section limits for the process and some model-dependent 
parameter space exclusion plots. 



Chapter 2 
Theory 



2.1 The Standard Model 

The experimental discoveries and theoretical developments that have led to our 
current understanding of the physical world are presently described within a single 
theoretical framework - the Standard Model of particle physics. It is the culmi- 
nation of a history of division and of unification, both driven by discovery. The 
division has been in the understanding of matter, as entities once thought funda- 
mental have successively been shown to be systems of even smaller entities. The 
atom (supposedly indivisible by definition) was found to be a system of negatively 
charged electrons with a positively charged nucleus. The nucleus was found to 
consist of positively charged 'protons' and neutral 'neutrons', and both these have 
been found to be comprised of yet smaller particles called quarks. At the same 
time understanding of the behaviour and interactions of these particles has been 
steadily consolidating. Maxwell showed that electricity and magnetism, previously 
thought to be entirely different phenomena, were in fact different manifestations of 
a single force subsequently dubbed electromagnetism. His theory also tied visible 
light and all other forms of electromagnetic radiation from radio waves to X-rays 
into a tight theoretical framework. Einstein provided the theories to link the worlds 
of the very fast and very large with the everyday world of Newton; and the other 
pioneers of quantum physics did the same for the very small. 

The principal discovery tools in particle physics this century have been particle 
accelerators/colliders. These have allowed the study of particle collisions at high 
energy, and unveiled a rich particle spectrum beyond that from which the everyday 
world is made up. These particles have been categorised by their properties. A 
principal distinction that can be drawn is based on particle spin. Spin is the 
'internal' angular momentum that a particle possesses, not by virtue of any motion 
it may have, but as an intrinsic property tied up in its very identity. All particles 
have spin in integer units of \ (in natural units). Fermions possess an odd multiple, 
bosons an even multiple. As a consequence, fermions and bosons have drastically 
different behaviour since fermions are subjected to the Pauli exclusion principle 
while bosons are not. Fermions are further split into leptons and quarks. The 
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distinction between these two groups is based on the forces they experience. 

There are four known forces which, in order of weakest to strongest, are gravity, 
the weak force, electromagnetism, and the strong force. All particles experience 
gravity, but since it is so weak at the microscopic level it has no observable effect on 
the interactions of individual particles. All fermions also experience the weak force. 
Particles will experience the electromagnetic force if they have electric charge. The 
distinction between leptons and quarks is that quarks experience the strong force 
whereas leptons do not. The 'charge' of the strong force is called colour. Particles 
can be red, green, blue or colourless. The nature of the strong force means that 
particles can only exist in isolation of others if they are colourless, i.e. if they 
contain equal amounts of red, green and blue quarks, a quark and an anti-quark 
of the same colour, or no quarks at all. Both quarks and leptons can be further 
split into 'up-type' and 'down-type'. All leptons/quarks of the same type have the 
same charge. They can also be split into three 'generations'. The first generation 
consists of the particles that make up everyday matter. The second and third are 
very similar to the first but heavier and unstable. The three down-type leptons 
are the electron, muon and tau. The three up-type are their associated neutrinos. 
The three down-type quarks are the down, strange and bottom, while the up-type 
are the up, charm and top. Each fermion has an associated anti-particle possessing 
the same mass and spin, but with opposite charge and colour (the anti-particle of 
the charge — | red up quark is the charge +| anti-red up anti-quark; note that 
red + anti-red = colourless). 

Since the earliest days of physics as a subject in its own right, it was under- 
stood that net amounts of certain quantities could neither be created nor destroyed 
- they are conserved. This was first observed in dynamical quantities such as en- 
ergy, momentum and angular momentum. But as research into matter continued, 
new conserved quantities were found that did not relate to the motion of parti- 
cles, but to 'internal' properties that appeared bound up in the particle's identity, 
such as electric charge. Initially these conservation laws were just that - laws in 
their own right which had to be inserted into the description of the world by hand. 
However, it was realised that conservation laws can actually be inescapable results 
of symmetries. That is, if the universe possesses a given symmetry, then any sys- 
tem should remain invariant under the respective transformation of the symmetry. 
Once this is inserted into the description of particles, a new prediction can appear 
which requires an associated quantity to be conserved. For example, invariance 
of a system under translation leads to momentum conservation, invariance under 
displacement in time leads to energy conservation, and invariance under a change 
in the phase of the wave-function leads to charge conservation. As opposed to a 
new conservation law which involves the introduction of a new arbitrary rule in 
need of further explanation, a new symmetry simply says that something is not 
important in our description of matter (it is a 'gauge') and so can be considered a 
simplification. 

This principle was extended further to great success. It was shown that con- 
sidering the phase of the wave-function of a charged particle to be invariant under 
local transformation necessitates the existence of the electromagnetic force. The 
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associated transformations form the group, U(l), and the associated eigenstate 
is called the B. Similarly, invariance under a change from up-type to down-type 
fermions necessitates the weak force. These transformations form the group SU(2), 
and the weak eigenstates are the W 1 , W 2 and W 3 . And invariance under colour 
rotation necessitates the strong force, whose associated group is SU(3) and whose 
associated eigenstates are the gluons (of which there are eight corresponding to dif- 
ferent permutations of colour-anticolour). Thus these three forces have been shown 
to be the results of symmetries, or 'gauge invariance'. [] 

Gauge invariance predicts the gauge quanta to be massless. This is fine in the 
case of the electromagnetic force since its range is known to be infinite. For a 
force to be transmitted over an arbitrarily long distance, the force mediators (the 
gauge quanta) have to be able to exist for an arbitrarily long period of time, and 
so by the Heisenberg uncertainty principle, AEAt < h, have to have an arbitrarily 
small amount of energy, and this can only be achieved if the quanta are massless. 
It is also fine for the strong force, because although its range is only of the order 
of 10~ 15 m, its quanta (the gluons) possess colour and so feel the force that they 
mediate, which can explain their restricted range. But while the quanta of the 
weak force can also interact with each other, the nature of their interaction cannot 
explain their very short range (~10~ 17 m). This implies that the quanta are heavy, 
and so could just be fixed by introducing masses for the W 1,2,3 by hand, but this 
destroys the gauge invariance and leads to unresolvable technical problems. The 
solution is provided by the Higgs mechanism. Higgs predicted the existence of 
a new spin-0 (or 'scalar') field which couples to the Standard Model particles. 
Its coupling to the SU(2) gauge quanta makes their symmetric states unstable, 
causing a 'spontaneous' breaking of the symmetry as it is realised in nature, and 
observed states in which the symmetry is not manifest. This implies the existence 
of new, massless bosons called Goldstone bosons. But rather than new physical 
states, these can be shown to manifest themselves as longitudinal degrees of freedom 
for the gauge bosons, giving the bosons mass in a manner consistent with gauge 
invariance. The coupling of the Higgs field to the fermions also accounts for their 
masses. A by-product of the Higgs mechanism is the (massive) Higgs boson, which 
has yet to be discovered despite the best efforts of experimental physicists. 

The mass eigenstates of the gauge bosons (the photon, 7, the Z and the W^) 
do not end up in a one-to-one correspondence with the weak eigenstates, but as 
mixtures of them - 

7 = cos 9 W B + sn\9wW 3 , 
Z = - sin 6 W B + cos 6 W W 3 , 
W ± = (W 1 =F iW 2 )/\/2 , 

where 9w is known as the weak mixing angle and is k,2%°. The masses of the 
weak gauge bosons, and the expected mass of the Higgs boson, are of the order 



1 Gravity is also a result of gauge invariance, specifically invariance under local coordinate 
transformation. But this is not yet incorporated in the Standard Model. 
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Generation 


+l=up-type 
— 1 =down-type 


Name Symbol Mass Charge 

(MeV) 


Leptons 


1 st 
-L 


-1 

+ 1 


electron e 0.511 —1 
e-neutrino v e < 3 x 10~ 6 


nnd 


-1 
+ 1 


muon jj, 106 —1 
/x-neutrino < 0.19 


ord 
O 


-1 
+ 1 


tau t 1780 -1 
r-neutrino v T < 18.2 


Quarks 




-1 
+ 1 


down d 3 — 9 — | 
up u 1 — 5 +| 


2^nd 


-1 
+ 1 


strange s 75 — 170 — | 
charm c 1150 - 1350 +§ 


3 rd 


-1 

+ 1 


bottom b 4000-4400 -| 
top t 174000 ± 5000 +| 



Table 2.1: Fermion summary table. Only upper limits exist on the masses of the 
neutrinos. Whether they have any mass at all is still under investigation. The 
masses of the quarks are not well defined because, due to the nature of the strong 
force, they cannot be isolated. All these values are taken from 0. 



Name 


Symbol 


Mass 


Charge 


Spin 






(GeV) 






Photon 


7 








1 


Z-boson 


Z 


91.2 





1 


W-boson 


W ± 


80.4 


±1 


1 


Gluon 


9 








1 


Higgs 


H 


> 114.1 








Graviton 


G 








2 



Table 2.2: Boson summary table. The upper limit on the Higgs mass is from the 
combined LEP search for the Standard Model Higgs boson ||. The graviton has 
not yet been directly observed. 
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of 100 GeV, and this is referred to as the 'weak scale'. The properties of all the 
particles of the Standard Model are summarised in Tables |2.1| and [2.2| . The charges 
are given as multiples of that on the positron. Also shown is the graviton, which is 
the hypothetical quantum of the gravitational force, which is not yet encompassed 
by the Standard Model and has not been directly observed. 



2.2 The hierarchy problem 

The Standard Model of particle physics has been an acclaimed success in describing 
with predictive power many aspects of the known particle spectrum. Despite this 
it is not regarded as the 'whole story'. It contains many free parameters which 
can only be determined by experiment, and it offers no explanation for many of 
the most intriguing aspects of the particle spectrum. Such as, why are there three 
generations? What are the origins of the values of the free parameters? Why is 
electric charge quantised in units of one third that on the electron? It certainly will 
not be able to describe physics at the Planck scale where quantum gravitational 
effects become important. And it is considered unlikely that in the many orders 
of magnitude between this scale and the weak scale there will be no intermediate 
new physics. Although questions such as these call into doubt the completeness of 
the Standard Model, they do not indicate any problems with the predictions that 
the Model does make. However, there is one area in which the Standard model is 
considered to have a flaw. This is in the Higgs sector. 

The Standard Model Higgs field is a complex scalar, H, with a potential given 

by 

V = i?\H\ 2 + \\H\* , 

where \i is the Higgs mass parameter and /i 2 < is required for electroweak 
symmetry breaking. Then the Vacuum Expectation Value (VEV) for the field 
is given by (H) = y^—/j, 2 /2X. It is known from electroweak measurements that 
(H) = 246 GeV. Thus if the Standard Model Higgs mechanism is responsible for 
electroweak symmetry breaking, /i 2 must be of the order of —(100 GeV) 2 . But /i 2 
is subject to loop corrections from every particle which couples to the Higgs field. 
These loop corrections are proportional to the square of the energy scale at which 
the loop integral is cut off, which would be the scale at which some new physics 
enters the picture to prevent further corrections. In general, this scale could be 
anywhere between the weak and Planck scale. To assume that it is low enough 
to prevent disastrously large corrections to fi 2 , but high enough to have had no 
measurable consequences given the energies of today's accelerators, is to ask it to 
lie in an extremely narrow region of the complete allowable scale. Without a model 
to explain why this should be so, it would seem like a very unnatural assumption. 

It is possible to remove this divergence by dimensional regularisation || , but one 
is still left with corrections proportional to the squares of the loop-particle masses. 
So the largest correction will be of the order of the (mass) 2 of the heaviest particle 
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that couples to the Higgs field, even if this is only indirectly via other interactions 
0]. This means that to preserve /i 2 at around — (lOOGeV) 2 , there would have to 
be no new particles that couple, directly or indirectly, to the Higgs field in a very 
large energy range up to and including the Planck scale. This is considered equally 
unnatural. 

One could assume that the various terms that contribute corrections to fi 2 
manage almost to cancel such as to leave the right amount. But for the reasons 
given above it would be expected that the largest of these terms would be very 
large indeed, very possibly of the order of the Planck mass. It would seem very 
unlikely that terms this large would cancel to leave just a few hundred GeV. For 
this to be the origin of the value of fi 2 would mean that the parameters of the 
theory would have to be fine-tuned to many significant figures. 

A physically justifiable model would be one in which the correct size for /x 2 
falls out naturally, without requiring any parameter conspiracies or a veto on new 
physics at higher scales. It is the failure of the Standard Model to provide this that 
is referred to as the "naturalness" , "fine-tuning" or "hierarchy" problem. 

2.3 Super symmetry as a solution 

Supersymmetry is a symmetry that relates fermions and bosons. Because of the 
relative minus sign between fermion and boson loop contributions to the Higgs 
mass, such a symmetry offers an immediate hope for a natural explanation of 
cancellations between loop-correction terms. Indeed, it is considered the primary 
attraction of supersymmetry (or 'SUSY') that, as an extension to the Standard 
Model and at least in its unbroken form, these cancellations are required to be 
exact. 

A SUSY transformation alters the spin of a state by one half unit, turning a 
bosonic state into a fermionic one and vice-versa: 

Q\boson) = | fermion); 
Q\ fermion) = \boson). 

The required operators, Q and (the hermitian conjugate of Q), must be anti- 
commuting spinors. Their algebra can be expressed in terms of commutation and 
anti-commutation relations; 

{Q,Qt} = P M ; 

{Q,Q} = {Q t ,Q t } = 0; 
[P»,Q] = [P»,Qi]=0. 

SUSY eigenstates which can be transformed into one another by the SUSY gen- 
erators are called "superpartners" . The collections of single particle states which 
form the irreducible representations of the SUSY algebra are called "supermulti- 
plets", and consist of pairs of superpartners. Since the (mass) 2 operator (— P 2 ) 
and the gauge group generators commute with Q and Q\ it follows that super- 
partner pairs have the same masses, gauge quantum numbers and couplings as 
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each other. It also follows (though not so simply) that each supermultiplet must 
contain an equal number of fermionic and bosonic degrees of freedom (ub = np). 
It is then guaranteed that all the loop contributions to the Higgs mass from a single 
supermultiplet must cancel exactly, and since under SUSY every particle is in a 
supermultiplet this forms a very neat solution to the hierarchy problem. 



2.3.1 SUSY and grand unification 

Supersymmetry also offers a solution to another problem of the Standard Model, 
that of grand unification, or lack thereof. 

Using the renormalisation group equations the Standard Model gauge couplings 
can be run up to higher scales, and it is found that they almost converge on a single 
value. This would be very desirable since it would be a clear indication of more 
fundamental physics at a higher scale. But measurements have shown that 'almost' 
is as good as it gets, they do not fully converge ||. The evolutions of the couplings 
depend on the particle content of the theory under consideration however, and 
so the additional particle content that results from supersymmetry changes the 
evolution. It just so happens that if no more supersymmetric particles than are 
absolutely necessary are included, the couplings evolve such as to unify at around 



10 GeV to within current experimental errors (see Figure |2.1| z |) . This could be a 



numerical coincidence, but it is certainly another motivation for supersymmetry. 

It should be stated at this point that there are many different models for the 
implementation of supersymmetry in nature, some more exotic than others. All 
require that the existing particle spectrum be at least doubled in size. Only the 
Minimal SuperSymmetric Model (or MSSM) is considered here. This is the model 
which introduces as few as possible new particles whilst still providing a consis- 
tent, phenomenologically viable model. As such it has been the most intensely 
investigated of all supersymmetric theories. 



2.4 The SUSY eigenstates 

How then, do the Standard Model particles fit into supermultiplets? Let's start 
with the fermions. Since n# = %, a supermultiplet containing a single massless 
fermion (np — 2) should also contain either a single massless spin-1 boson (% = 2) 
or two scalar bosons (2x[ne = 1])- In fact, for the theory to be renormalisable, 
the only spin-1 bosons allowed are the gauge bosons. Since no gauge boson shares 
the same quantum numbers as any fermion we are forced to propose a completely 
new set of scalar particles for the supersymmetric partners of the fermions. The 
two mass degenerate scalars associated with each fermion are the supersymmetric 
partners of the left-handed and right-handed states of the fermion. The naming 
convention for these new scalar particles is simply to append an 's' (for scalar) to the 
front of the name of their Standard Model counterpart. The notation convention 



2 Thanks to Aran Garcia-Bellido for providing me with this figure. 
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Figure 2.1: The evolution of the gauge coupling constants with energy, as predicted 
by both the Standard Model particle content alone (upper), and with the minimum 
additional content that could come from supersymmetry (lower). The energy scale, 
Q, is in units of GeV. Supersymmetry appears to allow grand unification. 



for all the new supersymmetric particles is that of a tilde above the letter used 
to denote the Standard Model equivalent. Thus we may refer for example to the 
"selectron" (e), "smuon" (p) and "sneutrino" (£>), or more generically to a "slepton" 
(I). Alternatively there is the "stop" (t) and the "sbottom" (6), both examples of a 
"squark" (q). These are all "sfermions" (fs). With the exception of the sneutrinos, 
there are two of each corresponding to the left-handed and right-handed states of 
the respective Standard Model particle. Thus the difference between a left-handed 
and a right-handed sfermion is nothing to do with spin states of the sfermions 
themselves (they are scalar), but instead refers to the spin-states of their Standard 
Model superpartners. 

Each gauge boson (n# = 2 before electroweak symmetry breaking) must be in 
a supermultiplet with a completely new massless spin-| fermion (again, a spin-| 
fermion would lead to an unrenormalisable theory). The naming convention for the 
fermionic partners of bosons is to add 'ino' to the end of the name of the Standard 
Model counterpart, replacing 'on' where it exists. Thus we have the "photino" (7), 
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"zino" (Z), "wino" (lU 1 * 1 ) and "gluino" (g), all of which are "gauginos". 

In supersymmetric theories it turns out that one Higgs boson is not enough. In 
order that all the Standard Model fermions gain mass, and in order to avoid triangle 
gauge anomalies, there must be at least two Higgs doublets ||. This doubles the 
number of real degrees of freedom to eight. After electroweak symmetry breaking, 
as in the Standard Model, three of these become the longitudinal modes of the 
Z and W^, although now five as opposed to just one are left over to form new 
scalar bosons. Thus supersymmetry dictates that there should be at least five 
Higgs bosons. Of these, three are neutral, one charge positive and one charge 
negative. Their spin-| superpartners are referred to as the "Higgsinos". All the 
supersymmetric partners of the Standard Model particles are collectively known as 
"sparticles" . 



2.5 R- parity 

Supersymmetry opens the possibility of Lagrangian terms which violate lepton 
and baryon number conservation. Although allowed they generally give rise to 
counterfactual predictions such as fast proton decay. Thus strong experimental 
limits exist on the couplings; they must either be very small or be forbidden by 
some conservation law. There are problems associated with the introduction of 
baryon and lepton number conservation as fundamental conservation laws in their 
own right. In the Standard Model they are satisfied 'accidentally' since there are no 
renormalisable terms which can violate them. In the MSSM they are not and their 
introduction as basic postulates would be a step backwards. Also they are known 
to be violated by non-perturbative electroweak effects. Instead, a new property 
known as "R-parity" can be introduced. This is defined for a (s)particle with 
baryon number, B, lepton number, L, and spin, S, as 

P R = (_ 1 )3(B-L)+2S _ 

It is +1 for all particles, and —1 for all sparticles, and there are no technical prob- 
lems associated with its introduction as a fundamental, multiplicatively conserved 
quantum number. Thus it can be assumed that Lagrangian terms are only allowed 
if the product of the R-parities of the involved fields is +1. This would rule out all 
baryon and lepton number violating terms. It would also have a huge impact on 
the phenomenology of SUSY. It would mean that sparticles can only be produced 
in pairs, and that all subsequent supersymmetric decay chains must end in the pro- 
duction of the lightest supersymmetric particle (LSP), which remains stable. Thus 
the LSP may be a good dark-matter candidate, and R-parity conserving models 
are subject to cosmological constraint. 

R-parity conserving and violating models of SUSY thus form two very con- 
trasting groups from the experimentalist's viewpoint. The search documented in 
this thesis assumes conservation of R-parity and so R-parity violation will not be 
discussed further. 
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2.6 Supersymmetry breaking 

As previously mentioned, supersymmetry dictates that sparticles should be degen- 
erate in mass with their corresponding particles. But if this were true they would 
have been discovered a long time ago. So if we are to continue with the hypothesis 
that supersymmetry is realised in nature, we must accept that it is a broken sym- 
metry resulting in sparticle masses significantly greater than their corresponding 
particles. However, the mass-splitting cannot be too large without unacceptably 
large loop corrections creeping back into the Higgs mass parameter. So if we are 
to retain supersymmetry as a solution to the hierarchy problem we find that 

IT^sparticle ^particle ^$ C?(l TeV) , 

leaving the sparticle spectrum still accessible to experiment. This is called 'soft' 
supersymmetry breaking, and it means we can write the effective Lagrangian of the 
MSSM as the Lagrangian of unbroken supersymmetry plus a Lagrangian containing 
the SUSY violating terms (the "soft" terms), 

£>MSSM — £>susy + £>soft ■ 

We do not need to make any assumptions about the origin of C so ft, our ignorance 
of the mechanism by which supersymmetry is broken can be parameterised by 
introducing every allowable soft term |7|]. This, however, is costly. While the 
unbroken MSSM introduces only one new, free, non-trivial parameter (the ratio of 
the VEVs of the two Higgs doublets, tan/5), the breaking of the MSSM in its most 
general form introduces 105 |§. Thus what was a highly predictive theory becomes 
so flexible it is almost untestable by experiment. It's clear then, that a model for 
the SUSY-breaking mechanism would be very desirable. 

If any random set of values were chosen for the 105 new soft parameters, one 
would almost certainly end up with a model incompatible with experimental obser- 
vation. Loop diagrams involving sparticles would lead to decays such as /i — > and 
b — > S7 with branching ratios far in excess of limits set by searches for flavour chang- 
ing neutral currents (FCNC's). And mixing between down and strange squarks 
would give extra contributions to K° — K° mixing incompatible with experimental 
measurements of the CP violating parameters. A desirable model would be one 
which naturally forbids all these counterfactual predictions. 

SUSY-breaking by renormalisable interactions at tree level does not do this. It 
also runs into technical problems when trying to give squarks and sleptons enough 
mass (to be heavy enough to have remained undetected thus far), and gauginos 
any mass at all. So it is generally assumed that the soft terms arise radiatively. 
That is, there exists a "hidden sector" of hitherto unhypothesised particles in which 
supersymmetry is broken. These do not couple to the "visible sector" of the MSSM 
particles and sparticles, but do share some interactions with it. Supersymmetry 
breaking is then communicated to the visible sector by loop corrections to the 
masses involving the gauge particles of the shared interactions. Additionally, if 
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these interactions are flavour-blind they can offer a natural explanation for the 
suppression of CP violation and FCNC's. 

2.6.1 The mass eigenstates 

The breaking of supersymmetry opens the possibility that the supersymmetry 
eigenstates do not correspond to the mass eigenstates, and therefore the possi- 
bility of mixing between sparticle states with the same quantum numbers. Thus in 
an analogous way to the mixing of the electroweak eigenstates following electroweak 
symmetry breaking, the gaugino and higgsino supersymmetry eigenstates can mix 
following supersymmetry breaking. The photino, zino, and neutral Higgsinos mix 
to form the "neutralinos" of which there are four, Xi 2 3 4- The wino and charged 
Higgsinos mix to form the "charginos" of which there are two, Xi2- The subscripts 
in each case indicate the mass order of the states, i.e. 

m^o > m^o > m^o > m^o . 

A4 A3 A2 Al 

In general sfermions with the same quantum numbers can also mix, but un- 
restricted mixing leads to the problems with FCNC's and CP violation already 
discussed. Thus any acceptable model for the origin of SUSY breaking must pre- 
dict that mixing between sfermions is kept slight. Most predict that it is negligible 
for the first two families, but possibly significant for the third. Thus the two 
mass eigenstates of the stau (and sbottom and stop) are differentiated by subscript 
numbers indicating which is the more massive, rather than by specifying a direct 
one-to-one correspondence with the Standard Model partner's spin-states. The re- 
lationships between the Standard Model particles, the SUSY eigenstates and the 
mass eigenstates are summarised in Table |2.3| . 

Of most importance in collider physics are of course the lightest of each type 
of sparticle. Thus, for the sake of brevity, terms like the neutralino, the selectron 
and the stau will be used in this thesis, where the is short for 'the lightest'. 

2.6.2 The gravitino 

The supersymmetry discussed thus far has been global supersymmetry. When we 
ask it to be a local symmetry, we find the requirement of invariance under local 
SUSY transformations implies invariance under local coordinate transformations, 
and this is the underlying assumption of general relativity. Thus local SUSY in- 
corporates gravity, and is called "supergravity" . 

The breaking of supersymmetry implies the existence of a spin-1/2 fermion, 
the "goldstino". Via the super-Higgs mechanism this is swallowed by the spin-| 
superpartner of the spin-2 graviton, the gravitino, giving it longitudinal degrees 
of freedom and thus mass. Without for the moment speculating as to the exact 
mechanism of SUSY-breaking, if we assume that its ultimate source is a VEV, F, 
in the hidden sector then by dimensional analysis we expect that the gravitino mass 
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Table 2.3: Summary of particle-sparticle correspondence 



will be of the order of 



2 



where M R = 2.4 x 10 18 GeV/c 2 is the reduced Planck mass. This is because the 
gravitino must lose its mass in the limits of vanishing SUSY-breaking (F — > 0) 
and/or vanishing gravitation (Mp — > 00). The exact relationship is 



ni3 



There are two candidates for the interactions that mediate SUSY breaking to 
the MSSM - gravity and the ordinary gauge interactions. The phenomenological 
importance of the dependence of the gravitino mass upon F will become clear in 
the following discussion of these two mechanisms. 



2.6.3 Gravity mediated supersymmetry breaking 

The historically favoured mechanism for the breaking of supersymmetry is via 
gravity 00. In this scenario SUSY-breaking is transmitted to the visible sector 
by gravitational interactions, or more generally, by interactions of gravitational 
strength. Then by dimensional analysis the soft terms should be of the order of 

F 
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suggesting that, for m so f t of the order of 100 GeV, vF should be of the order of 10 10 - 
10 11 GeV, giving a gravitino mass comparable to m so ft- Thus the LSP is almost 
certainly not the gravitino. Since, under R-parity conservation, the universe will 
have a relic density of LSP's left over from the Big Bang, the LSP almost certainly 
cannot be charged without having already been discovered or producing a universe 
incompatible with current observations. Thus the LSP is most often assumed to be 
the neutralino, and sometimes the sneutrino. The neutralino is favoured because 
most models predict it to be the LSP, and because there are problems with creating 
sneutrino-LSP models which predict a relic LSP density high enough for it to form 
a realistic dark- matter candidate |nj . This is because sneutrinos annihilate easily 
via the Z since their coupling to it is full strength, whereas the lightest neutralino, 
which most models predict to be mainly bino (the supersymmetric partner of the 
B, see Page O), does not. 



2.6.4 Gauge Mediated Supersymmetry Breaking (GMSB) 



Under gauge- mediated supersymmetry breaking ( pT]| and, for reviews, p2| [^3|| ) it is 
the normal gauge interactions that communicate the breaking to the visible sector. 
The breaking is communicated from the hidden high-energy sector down to the 
visible sector by a set of chiral supermultiplets which couple to the breaking sector 
and also indirectly through gauge interactions to the visible sector of the MSSM. 
Now if the mass scale of supersymmetry breaking is denoted F, and the mass scale 
of the messenger particles is denoted M mess , then by dimensional analysis 



m 



soft 



and so for m so f t of the correct order of magnitude a/F can be very much lower; 
as low as ~ 10 4 GeV since there is no reason why M mcss cannot be as low as 
possible without contradicting any current observations. Thus the gravitino is 
almost certainly the LSP under GMSB and all supersymmetric decay chains will 
end in the production of a gravitino (which forms a good dark matter candidate). 
Simple GMSB models can be described by six parameters : 



F - the scale of supersymmetry breaking in the messenger sector 

M mess - the messenger mass scale 

N5 - the number of complete SU(5) multiplets that make up the mes- 
senger sector 

A - the universal mass scale of the SUSY particles 

tan p - the ratio of the VEVs of the two Higgs doublets 

sign(/i) - (fji/\[ji\) where \i is the higgsino mass parameter 



The messenger sector is assumed to consist of a number of complete SU(5) mul- 
tiplets because then it will not interfere with the running of the coupling constants 
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to the extent that the MSSM prediction of grand-unification is scuppered. But 
unification can still be lost for a sufficiently large N5, so it is normally assumed to 
be no more than around five. The most important of these parameters is A, since it 
is principally responsible for the magnitude of all the gaugino, slepton and squark 
masses. 

The gaugino mass parameters at the messenger scale are given from the one-loop 
diagrams as 

m a = ^-N 5 A , 

4:71 

and the scalar masses from the two-loop diagrams (they do not receive corrections 
at one-loop order) as 

where C3 = | for colour triplets, for colour singlets; C2 = § for weak doublets 
and zero for weak singlets; and Y is hypercharge (Q = T3 + \Y). Thus the hier- 
archy, a 3 > a 2 > ai, means that strongly interacting sparticles are more massive 
than weakly interacting ones. Since the breaking mechanism is flavour blind, all 
sparticles with the same electroweak quantum numbers are automatically degen- 
erate in mass. Radiative corrections can in general cause mixing between the left- 
and right-handed sfermions, but the amount of mixing depends on the respective 
Yukawa couplings and so they can be neglected for at least the first two genera- 
tions. Thus gauge-mediation offers a natural suppression of FCNC's. Mixing in 
the third generation increases with tan (3 and causes the lightest stau to become 
lighter than the selectron and smuon. In evolving these masses down to the weak 
scale using the renormalisation group equations, the masses also gain a logarithmic 
dependence on M mess . 

A principal difference between models then, and of great importance to the 
manifestation of SUSY in detectors, will not be the identity of the LSP, but of 
the NLSP (the next-to-lightest supersymmetric particle). This must be either the 
neutralino, the right-handed sleptons collectively, the lightest stau in the case of 
large mixing, or a combination of these. The analysis described in this thesis 
assumes either the stau, or the three sleptons collectively. 

Since the gravitino inherits the non-gravitational interactions of the goldstino 
it absorbs, decay rates to it can be high enough for it to play a role in collider 
physics. The lifetime of a slepton NLSP, tj, depends on the gravitino mass and is 



m 2 sca iar = 2iV 5 A 2 
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So depending on the gravitino and slepton masses the slepton lifetime can be such 
that the decay-length is anywhere between so short it is unobservable in a detector 
to so long it will escape a detector before decay. An almost identical equation holds 
for the neutralino if it is the NLSP, and so the same applies. 



2.6.5 GMSB search signatures at LEP 

This section outlines the various ways that GMSB could be detected at LEP. It 



is basically a summary of [13], in which the phenomenology of a general class of 
GMSB models is considered. 

Search strategies must be geared towards hypotheses for the identity of the 
NLSP. Because the gaugino, squark and slepton masses are roughly proportional 
to their gauge charges only the lightest neutralino and sleptons form plausible 
NLSP candidates. Since the gaugino masses are proportional to N5, while the 
squark and slepton masses go only as \/N5, the neutralino and the sleptons are 
favoured as NLSP's for low and high N5 respectively. 

In the neutralino-NLSP scenario, GMSB could be discovered through the direct 
production of neutralinos and their subsequent decay to a photon plus gravitino. 
Neutralinos can be pair-produced either through e + e~ annihilation via the Z, or 
via t-channel selectron exchange (see Figure ]2.2|) . The t-channel is generally the 
more dominant because of the relatively large coupling. The detector signature 
of such events is a question of the neutralino decay-length. In the case of very 
short decay-length, the signature is two acoplanar photons (i.e. the two photon 
directions and the beam axis will not lie in a common plane) and missing energy 
(i.e. the energy of the event will be less than the LEP centre of mass energy). The 
missing energy is that carried away by the gravitinos - being neutral and colourless 
they will be invisible to the detectors. If the decay-length is of the order of a typical 
detector size (i.e. ~lm) then the signature is still acoplanar photons with missing 
energy, but now it can be possible to detect the fact that the photons are not 
originating from the primary interaction point (the IP) which will greatly enhance 
the discovery potential because backgrounds (which should involve photons from 
the IP) will be easier to remove. Should the decay-length be too long however, 
the neutralinos will escape the detector before decay leaving no evidence that any 
event actually occurred, becoming an invisible channel with no discovery potential 
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at LEP. 




Figure 2.2: Neutralino production and decay topologies. Note that 'x' refers to Xi 
- an abbreviation that will hold throughout the rest of this document. The two 
diagrams on the left show production in the s and t channels. The two on the right 
show decay in the neutralino-NLSP scenario (upper) and slepton-NLSP scenario 
(lower). 

If one or all sleptons are light enough, they can contribute to neutralino pro- 
duction under the neutralino-NLSP scenario through cascade decay. Sleptons can 
be pair-produced through annihilation via the 7 or Z, and selectrons only can also 
be produced through the t-channel exchange of a neutralino (see Figure |2.3| ), al- 
though the s and t channels generally interfere destructively in this case. They can 
then each decay to their respective lepton plus a neutralino (which will be favoured 
over decay directly to a gravitino in the neutralino-NLSP scenario). The signature 
here then will be two acoplanar photons and missing energy, plus two (typically 
low momentum, or 'soft') leptons from the IP. However the slepton cross-section 
under the neutralino-NLSP scenario is expected to be much smaller than that of 
the neutralino, and so this is not of great importance. 




Figure 2.3: Slepton production and decay topologies. The two diagrams on the left 
show production in the s and t channels (only selectrons can be produced in the 
t-channel). The two on the right show decay in the slepton-NLSP scenario (upper) 
and neutralino-NLSP scenario (lower). 

There are two slepton-NLSP scenarios, depending on the amount of stau mixing. 
In the case of little mixing the selectron, smuon and stau are degenerate in mass 
(and so also in lifetime). In this case they act collectively as the NLSP, and this 
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is called the slepton co-NLSP scenario. For large mixing (high tan/5), the lightest 
stau becomes significantly lighter than the selectron and smuon and can act as the 
sole NLSP. This is called the stau-NLSP scenario. Direct slepton pair-production is 
then going to be the most obvious process to look for in both these cases, with each 
slepton decaying to its respective lepton and a gravitino. If the decay-length is short 
this will happen at the IP and so the signature will be two acoplanar leptons plus 
missing energy (note that the same signature is important in searches for gravity 
mediated SUSY where the sleptons decay instead to stable neutralinos) . Again, if 
the decay-length is longer, the fact that the leptons are not originating from the IP 
can be detected, enhancing the signature in the same way as for the neutralinos. 
But now since the slepton is charged its track through the detector can be visible, 
and so therefore can the slepton decay vertex, enhancing the distinctiveness of 
the signal even further. If the decay-length is so long that the slepton escapes 
the detector before decay, then the sleptons will be the only particles seen in the 
detector and the full centre of mass energy will be observed. These events will look 
similar to e + e _ — > events, but a detector with sufficient single-particle mass 

resolution (through calculation of a/ E 2 — p 2 and/or measurement of the particle's 
rate of ionisation) will be able to distinguish them. 

If the neutralino is light enough, it can contribute to slepton production in 
the same way that the reverse can be true in the neutralino-NLSP scenario. But 
now the contribution through cascade decay in this way can be comparable to, 
and often greater than, direct pair production. Thus in the slepton (co or stau) 
NLSP scenario, slepton production though the cascade decay of neutralinos forms 
a good discovery channel. In this case the signature is similar to that from direct 
slepton production, but now with two additional softer leptons from the IP. Also, 
the two harder leptons need not have the same charge now since the decays of each 
neutralino will be independent. 

It is also possible that the stau, or all three sleptons, are sufficiently degenerate 
with the neutralino such as to form a slepton-neutralino co-NLSP scenario. In 
this case events from both neutralino and slepton pair production will contribute 
simultaneously, and cascade decays will be highly suppressed. 



2.6.6 Search statuses and limits 

This section briefly summarises the searches that had been performed, and mass 
limits that had been obtained, after the LEP data taking period of 1998 during 
which approximately 170 pb -1 of data were taken by each experiment at a centre 
of mass energy of 189 GeV. This marks what was essentially the starting point for 



the analysis described in this thesis. The current limits will be given in Section [77?. 

ALEPH performed searches for neutralino pair production under the neutralino- 
NLSP scenario for both (effectively) zero and observable neutralino decay-length 



1511 [M ; and for slepton pair production under the slepton co-NLSP and stau-NLSP 



scenarios for zero, observable, and very large decay- lengths |PT|| ||18||[P1?|| ; and for the 
cascade production of sleptons from the pair production of neutralinos under the 
slepton co-NLSP and stau-NLSP scenarios for zero decay-length pi. All these 
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searches were brought up to date with the 189 GeV data and had their results 



collectively interpreted in [16 



In the neutralino-NLSP scenario a lower neutralino mass limit of 91 GeV was 
obtained for zero decay-length, falling to 55 GeV for ct x < 100 m. In the stau-NLSP 
scenario a lower stau mass limit of 67 GeV was obtained for any lifetime. In the 
slepton co-NLSP scenario a lower limit on the common slepton mass of 84 GeV was 
obtained for any lifetime (under the assumption that selectron production proceeds 
only via the s-channel). In the case that the neutralino mass is less than 87 GeV 
and the neutralino-stau mass difference is greater than the tau mass and that the 
stau lifetime is negligible, the search for slepton production by cascade decay from 
neutralinos allowed the limit on the stau mass in the stau-NLSP scenario to be 
increased to 84 GeV (under the assumption that the neutralino is mainly bino and 
that its mass is | that of the selectron). 

The results from all these searches were interpreted in terms of the (\/F, M mess , 
N5, A, tan/3, sign(/i)) parameter space according to the GMSB model described 
in [20|. A scan over the parameter space used the results to determine the regions 
excluded. It provided a lower limit of 45 GeV on the NLSP mass for any NLSP 
lifetime under any scenario, a lower limit of 9 TeV on A, and of 2 x 10 -2 eV on the 
gravitino mass. 

Similar work was performed by DELPHI and OPAL in [^1|] and respectively, 
but the former does not consider the neutralino-NLSP scenario, and the latter 
deals only with zero NLSP decay-length (although a slepton-with-lifetime search is 
detailed in IE3V 



Chapter 3 

The ALEPH detector at LEP 



3.1 The LEP collider 



'LEP' is the Large Electron-Positron collider at CERN. It is a circular device of 
diameter 8,486 m lying in a tunnel under the Swiss- French border near Geneva. 
During operation, electrons and positrons are accelerated through an evacuated 
beam-pipe in opposite directions by RF cavities. Their roughly circular orbits 
within the machine are created by bending the particle beams with dipole mag- 
nets. This inevitably leads to energy loss through synchrotron radiation, which is 
replaced by the cavities. This is the reason for the large size of LEP as the amount 
of energy lost in this way is inversely proportional to the accelerator's radius. The 
counter-propagating particles are allowed to collide with equal and opposite mo- 
menta at four points, each surrounded by a detector whose purpose is to observe 
and record new particles produced by the collisions. Each detector constitutes 
a separate experiment, these are ALEPH, OPAL, DELPHI and L3. For general 



papers on these detectors see p4|,Pq|,[p6| and p7| respectively. 

The volume at the centre of each detector in which the beams overlap is known 
as the 'beam-spot' or 'luminous region', and it is within this that collisions occur 
and from where (primary) particles originate. Since the beams consist of discrete 
bunches of particles, the size of the luminous region is defined by the size of the 
bunches. These span 200 /iin horizontally and 8 /im vertically in the plane perpen- 
dicular to the beam direction, and 1 cm parallel to the beam direction. 

In the first stage of its operation, known as LEP1, LEP ran such that the 
centre-of-mass energy of its colliding particles was equal to the mass of the Z boson 
(~ 91 GeV), allowing it to fulfil its first goal of performing detailed studies of the Z. 
Subsequently it was upgraded and its energy increased through LEP1.5 to LEP2, 
passing the threshold for, and allowing the study of, W boson pair production and 
reaching in excess of 200 GeV. 
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3.2 An overview of ALEPH 



ALEPH is a general purpose detector which allows the study of almost all types 
of e + e~ interactions. It surrounds the interaction point with near An solid angle 
coverage, as shown in the cut-away diagram in Figure |3. 1| . Its main bulk consists of 
several different subdetectors forming coaxial jackets about a 5.3 cm-radius beam- 
pipe, centred on the interaction point (IP). 
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Figure 3.1: The ALEPH detector. 

Particles produced at the IP and travelling outwards through the detector first 
encounter a series of three, low-density tracking subdetectors designed primarily 
to give information on the trajectories of charged particles. Since the tracking 
volume lies within a strong magnetic field, the curvature of a charged particle's path 
through the field yields its momentum. They then reach a high-density calorimeter 
arrangement designed to bring them to rest and, in doing so, make a measurement 
of their energy. Any particles that manage to penetrate this, pass finally through 
the muon chambers designed to tag muons leaving the detector. The magnetic 
field is produced by a superconducting coil located between the electromagnetic 
and hadron calorimeters, and runs parallel to the beam-axis with a strength of 
1.5T. 

The Cartesian (x, y, z) coordinate system used to describe positions within 
ALEPH is defined as follows. The origin is the geometrical centre of the detec- 
tor, also the nominal IP. The z-axis runs parallel to the beam-pipe in the e~ beam 
direction. The x-axis points towards the centre of LEP (to the right as seen by 
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an incoming positron), while the y-axis points vertically upwards^. A cylindrical 
(r, <f), z) coordinate system is also used, in which the origin and z-axis coincide with 
those of the Cartesian system. The planes = 0° and <p = 90° contain the x-axis 
and the y-axis respectively, and r = \/ x 2 + y 2 . In addition, it is often useful to 
refer to the angle made between the line that joins a point to the origin and the 
z-axis, and this is denoted by 9. 

A detailed description of the performance of the ALEPH detector is available in 



28] . It should be noted that ALEPH has been undergoing constant modification 
over the years it has been operating, from minor changes too slight to be of note, 
to major changes like the replacement of the vertex detector in 1995. This analysis 
however, is only concerned with the data-taking period, 1997 onwards. Only the 
state of the detector in this period is described here, and when this state differs 
from previous years this will not necessarily be noted. With this in mind, there 
now follows a detailed description of the main components of ALEPH. 



3.3 The Vertex DETector (VDET) 

The VDET is the innermost subdetector of ALEPH. It is a silicon-strip detector 
which provides high-accuracy positional information on the paths of charged parti- 
cles close to the IP both in the rep plane and the z direction. This allows accurate 
determination of their trajectories in this region and thus gives information on the 
point of origin of a charged particle, whether that be from the IP or not. This is of 
course most useful for the study of B-physics, allowing the accurate reconstruction 
of the decay vertices of B-mesons close to the IP. 

The active part of the VDET is made up of 24 'faces' arranged in two, approx- 
imately cylindrical layers. The inner layer consists of 9 faces and is at the lowest 
possible radius allowed by the presence of the beam-pipe (~ 6.3 cm), the outer 
consists of the remaining 15 and is at the greatest possible radius allowed by the 
presence of the Inner Tracking Chamber (~ 10.3 cm) - thus maximising the lever 
arm. The faces overlap by ~ 0.2 mm in rip to ensure complete coverage. 

The active units of the VDET are the 'wafers', these are double-sided silicon 
strip detectors of size 52.6 mm x 65.4 mm x 0.3 mm. Essentially junction diodes 
operated in reverse bias, each is an n-type substrate with 1021 p+ readout strips on 
the junction {rip) side and 640 orthogonal readout n+ strips on the ohmic (z) side. 
The passage of a charged particle through the wafer generates ionisation which 
is picked up as signal by nearby strips. Clusters of signals are combined using a 
'centre-of-gravity' algorithm to produce VDET hits. 

Three wafers and a support for the readout circuits all glued together form a 
'module', two modules mounted with their readouts at opposite ends on two beams 
form a face, active length 400 mm, total length 500 mm. One beam is made of 
Kevlar epoxy and serves to electrically insulate the z side of the modules, the other 



^his is only approximate. The slight tilt of the LEP ring with respect to the horizontal means 
that there is a small angle between the j/-axis and the vertical. 
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is made of Carbon fibre epoxy and provides the necessary mechanical strength of 
the face. The faces are then mounted at each end on carbon fibre support flanges, 
which are joined together by a carbon fibre cylinder which sits between the two 
active layers and consists of two thin skins spaced 20 mm apart by a corrugated 
web. The whole support structure is water cooled. The readout cables and water 
pipes are supported by a second pair of 20 fim thick aluminium flanges, joined by 
two 200 /im thick carbon fibre cylinders which serve both to support the flanges and 
as inner and outer protective skins for the VDET. Figure |3]^ shows a cross-section 
of the completed structure in the rz plane and the face geometry as seen in the r0 
plane, Figure |3]3| shows details of a face as seen from above and below and a three 
dimensional view of the fully mounted VDET. 




Figure 3.2: Cross-sections of the VDET showing structural composition in the rz 
plane (left) and face geometry in the r</> plane (right). 




Figure 3.3: A VDET face (left) and the fully mounted VDET (right). 
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A charged particle passing through a face well inside its active region stands a 
99% chance of leaving a reconstructed hit in both the r<j) and z views. The chance 
of that hit being assigned to its reconstructed track in one view when it has been 
assigned in the other, giving a three-dimensional hit for the track, is 93% for the 
z view and 95% for the r<p view. The spatial hit resolution for primary tracks at 
9 = 90° is 10 /mi in the rip view, with no significant change for decreasing 9, and 
15 /mi in the z view, worsening to 50 /mi by |cos(#)| = 0.85. The 40 cm active 
length of the VDET (virtually) guarantees at least one hit up to | cos(0)| = 0.95. 

3.4 The Inner Tracking Chamber (ITC) 

The ITC is a cylindrical wire drift-chamber providing tracking information in the 
region 26 cm > r > 16 cm and \z\ < 100 cm. It has a fast readout allowing it to be 
used in the Level- 1 trigger decision and provides the only tracking information used 
in that decision. For a charged particle passing through its full active thickness it 
can provide eight hits with a resolution of ~ 150 in r<p but only ~ 5 cm in z. 

The wires are strung parallel to the z-axis between two 25 mm thick aluminium 
end-plates. The end-plates themselves are connected by two carbon fibre tubes 
which form the cylindrical surfaces of the ITC. The outer tube bears the tension 
of the wires, and has an inner radius 57 cm, a thickness of 2 mm, and has a 25 /mi 
layer of aluminium foil on both surfaces to screen the chamber from RF interference 
and improve the uniformity of its field. The inner tube, which provides support to 
the VDET, is 600 /mi thick with a 50 /mi thick layer of aluminium foil on its outer 
surface for the same reason. The high voltage is supplied by distribution boxes 
mounted on electrical end-flanges 20 cm beyond the end-plates. These end-flanges 
plus the cylinders form a hermetically sealed volume which contains a gas mixture 
of four parts argon to one part carbon dioxide at atmospheric pressure. 

There are four different types of wire. The sense wires, responsible for attracting 
the ionisation produced by passing charged particles, are held at a positive potential 
of ~ 2 kV. Each is surrounded by five field wires and a calibration wire all held at 
ground potential, forming a hexagonal 'drift cell' (see Figure |3.4|) . Between each 
pair of drift cell layers is a layer of guard wires, these support circular hoops of 
aluminium wire whose purpose is to limit the damage that would be caused should 
any wire break. There are 96 sense wires per layer in the first four layers, and 144 
per layer in the outer four layers. 

An ITC hit is created when a charged particle passes through a drift cell. Its 
passage ionises the gas in its immediate vicinity causing a pulse of negative charge 
to drift to the sense wire. The subsequent current pulse is detected at both ends. 
The distance of the particle from the wire is obtained by converting the drift time 
into a drift distance using a parameterisation of the non-linear relationship between 
the two. ITC hits have an inherent r<p ambiguity due to the azimuthal symmetry 
of the drift cells, the result of which is that each hit can represent one of two 
points on either side of the relevant sense wire. This ambiguity can be removed in 
practice since drift cells of neighbouring layers are offset from each other by half 
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Figure 3.4: A detail of an ITC end-plate (left) and the drift cell geometry (right). 

a cell width, thus only the correct set of possible hit points will line up to form a 
valid track. The z coordinate of a hit is obtained from the difference in the time 
of arrival of the pulse at each end of the sense wire. Note that no more than one 
hit can be assigned to each wire in a given event. 

3.5 The Time Projection Chamber (TPC) 

The TPC is ALEPH's main tracking chamber. It accurately measures the tra- 
jectories of charged particles by providing three-dimensional hits at 21 separate 
radii in the range 170cm > r > 40cm and \z\ < 220cm. It also measures the 
spatial rate of energy loss due to ionisation (dE/dx) of a particle, which enhances 
particle identification by complementing information from the calorimeters. Struc- 
turally, the TPC consists of three main elements - the field cage (two cylinders, 
inner and outer), two circular end-plates and eight 'feet' (four attached to each 
end-plate, they transmit the weight of the whole structure, plus that of the ITC, 
to the magnet cryostat). A diagram is shown in Figure |3.5| . 

The cylindrical field cage is aligned with the z-axis and has an inner radius of 
31 cm and an outer radius of 180 cm. This, together with the end-plates, forms the 
gas-tight volume of the TPC which contains a mixture of argon (91%) and methane 
(9%) held at slightly above atmospheric pressure. The volume is vertically bisected 
by a circular mylar membrane coated in conducting graphite paint and held at a 
high negative voltage (typically — 27 kV). The end-plates are held near ground while 
electrodes along the inner and outer field cages are held at potentials such that the 
resulting electric field (of ~115 Vcm -1 ) is uniform and aligned with the z-axis. 

A charged particle passing through the volume ionises the gas along its path, 
and the electric field causes the resulting electron cloud - an image of the particle's 
trajectory - to drift towards the nearest end-plate, while its lateral diffusion is 
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Figure 3.5: TPC overall view. 



limited by the parallel magnetic field. 

The end-plates are divided into 'sectors' which are the active detecting regions. 
Each is a proportional wire chamber over a series of concentric 'pad' rows. The 
pads, 6.2 mm x 30 mm (5i(f) x 5r), provide the three-dimensional hit coordinates, 
while the wire chambers provide the dE/dx information. A diagram of an end-plate 
showing the sector, wire and pad geometry can be seen in Figure |3.6j . There are 
three different sector types (M, W and K). Significant dead regions, 24 mm wide, 
exist between neighbouring sectors' radial boundaries such that any portion of a 
particle's path in r0 that lies over a dead region will not produce hits. As such, the 
relative geometry of the three sector types was chosen such that the dead regions 
zigzag, limiting the number of hits a particle can lose in this way. 

Ionisation from a particle's path arriving at the end-plate encounters first a 
layer of wires called the gating grid (see Figure |3.7|) . If the Level-1 trigger has 
reached a 'yes' decision this will be held at a potential such that it is transparent 
to the passage of ionisationP]. Next it will encounter the cathode grid which is 
held at null potential to provide a shield between the uniform field of the main 
TPC volume and the field produced by the third and final layer. This is the sense 
grid and consists of alternating sense and field wires, the sense wires being held 
at null potential while the field wires are held at a high positive potential. As a 
result, incoming electrons avalanche towards the sense wires. The time of arrival 



2 i.e. the gate will be open. If there has been no such decision a potential difference of 200V will 
exist between neighbouring wires making the grid opaque, (gate closed). This is to stop positive 
ions generated in the avalanches described later from drifting back into the main TPC volume 
and causing distortions in the electric field. A Level-1 'yes' decision holds the gate open for 45 /is, 
the maximum drift-time of the TPC. 
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Figure 3.6: Overall geometry of a TPC end-plate showing sectors, wires and the 
21 concentric pad rows. 



and profiles of the resulting charge pulses are measured both by the sense wires 
and, through capacitive coupling, by the nearby pads. 

The integrated charge of each wire pulse is proportional to the dE/dx of the 
particle multiplied by the length of the particle's path that projects onto the wire. 
They are recorded and are used to calculate the dE/dx of each track once recon- 
struction has taken place. Clusters of pad pulses form the (r,(f>,z) TPC hits. Pad 
pulses which have multiple peaks in their time-profile are split, allowing more than 
one pulse per pad and thus ultimately, hits that overlap in r<p. The radius of a hit 
comes simply from the radius of the relevant pad row (clusters do not transcend 
rows). The value comes from a combination of the 0's of the relevant pulses, 
although the exact nature of the combination depends on the number of pulses 
involved. The z value is calculated from the cluster time using the drift velocity, 
the cluster time being a charge-weighted average of the time-of-arrival estimates of 
the pulses. The drift velocity is known from laser calibration. 

A single-hit resolution of 173 /im in r<p and 740 /im in z has been measured 
from leptonic Z° decays (which involve few, high momentum tracks leading to 
clear, accurately reconstructed events). Average resolution decreases with parti- 
cle momentum and dip angle (as such, two-photon events have about the worst 
resolution). 

The fractional momentum resolution, which is proportional to momentum, was 
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10mm 



Figure 3.7: Schematic diagram of a sector edge showing wires and pads. The 
strips on the left are held at potentials such as to reduce field distortions at sector 
boundaries. 

For tracks reconstructed purely from 
3 GeV _1 . 

3 GeV-\ 
and with VDET hits, to 

^ = 0.6 x 10- 3 GeV _1 . 

3.6 Material density within the tracking volume 

Clearly, from the point of creation of a particle to its entrance into the calorimeters, 
any interaction with the material of the detector is undesirable since it will alter 
the particle's 4-vector and so corrupt its measurement. Such interactions can also 
create new particles which may be confused with those resulting from the physics 
at the interaction point. This would be especially problematic close to the beam 
axis where information about a particle's momentum, or even existence before an 
interaction would be minimal or non-existent. An ideal tracking device would 



measured from ideal Z — > ji + fx eventsfj. 
TPC hits it was found to be 

Ap 

-5- = 1.2 x 1(T 

With the benefit of ITC hits this drops to 

Ap 

4 = 0.8 x 10" 



3 'Ideal' means the angle between the two tracks is greater than 179.7° and the total ECAL 
energy unassociated with the muons is less than 100 MeV. 
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therefore be transparent to the particles passing through it, i.e. of zero density. 
Obviously this is not achievable in the real world where the amount, distribution 
and density of material comprising a subdetector can only be minimised within the 
constraints of the design criteria of that subdetector, both as a measuring device 
of sufficient accuracy and as a physical structure of sufficient strength. 

The principal high-density components that exist within the ALEPH tracking 
volume are summarised in Table 3.1 , and it is these that are responsible for the 
vast majority of material interactions at non-negligible polar angles. Most of this 
material is necessary to isolate the various gas systems that exist within ALEPH, 
i.e. the virtual vacuum of the beam-pipe, the argon-carbon dioxide mixture of the 
ITC and the argon-methane mixture of the TPC. The total number of radiation 
lengths of material to be penetrated by a straight particle from the interaction point 
before reaching the first VDET layer and first TPC pad are shown in Figures |3.8| 
and respectively as a function of polar angle. The former is important since 
the first VDET layer is the first active part of the detector that a particle will 
typically encounter. Any interaction that happens at a lower radius happens 'in 
the dark', corrupting the relevant particle's 4- vector before it has a chance even to 
be detected. Also, any new particles that are created as a result of the interaction 
will be more difficult to reject since they will not be lacking any hits. The amount 
of material before the first TPC pad is of interest since there is no more solid 
matter between here and the outer edge of the tracking volume. Figure then 
serves as a summary of the problematic material listed in Table [3.1| as a function 
of polar angle, and shows that a high-energy photon emitted from the IP with 
9 = 90° stands a 6% chance of pair-converting before entering the TPC. Note that 
in general, as the polar angle decreases, the thickness of material to be penetrated 
increases, due to the decreasing angle of incidence into cylindrical components of 
the detector. Peaks and lumps are due to circular end-plates and non-cylindrical 
support structures. 



3.7 The Electromagnetic CALorimeter (ECAL) 

The 'ECAL' is a lead and wire chamber sampling calorimeter with a thickness of 22 
radiation lengths. It performs a multi-stage energy measurement of electromagnetic 
particles. Opaque to electrons and photons, it can yield their total energy. Hadrons 
and muons typically penetrate its full thickness. In measuring both the transverse 
and longitudinal shower profiles created by particles it also provides information 
on particle identity. 

Again, it is comprised of a barrel plus two endcaps aligned with the z-axis. 
Each is segmented into twelve modules in 0. A module (there is little difference 
whether it be in the barrel or an endcap) is a hermetically sealed unit containing 
45 'layers' in an 80%-xenon to 20%-carbon dioxide gas mixture held at slightly 
above atmospheric pressure. A layer consists of a lead sheet, an aluminium extru- 
sion containing an anode wire plane, and a plane of cathode pads separated and 



insulated from the wires by a graphite-coated mylar layer (see Figure |3.10| ). An 
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Figure 3.8: Material before the first VDET layer as a function of polar angle, 9. 
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Figure 3.9: Material before the first pad of the TPC as a function of polar angle, 9. 
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Component 


Purpose 


Material 

(or principal 
materials if many) 


Geometry 


Thickness 
(mm) 


Radiation 

lengths 
(7o ot X J 


Beam-pipe 


Holds LEP vacuum & 
supports VDET 


Be 


Cylinder, 
r = 5.3 cm 


1.1 


0.3 


VDET 


Active subdetector 
+ support structure 


Si, 

carbon fibre 
& Araldite 


("YimplpY ape 

sectKa 


N/A 


1.5 


ITC inner wall 


Contains ITC gas & 
supports VDET 


Carbon 
fibre 


Cylinder, 
r = 12.8 cm 


0.6 


0.3 


ITC end-plates 


Contain ITC gas & 
transmit wire tension 
to outer wall 


Al 


Holed discs, 
30 cm > r > 12.8 cm, 
z = ±100 cm 


25 


28 


ITC outer wall 


Contains ITC gas & 
counters wire tension 


Carbon 
fibre 


Cylinder, 
r = 28.5 cm 


2 


1 


TPC inner wall 


Contains TPC gas 

& forms inner 
electric field cage 


Al, mylar 
& Nomex 
honeycomb 


Cylinder, 
r = 31 cm 


10.5 


2.3 



Table 3.1: A summary of the principal high- density components within the tracking 
volume. Numbers quoted in the geometry column give the inner surface of the 
component. 



incoming photon or charged particle will generate an electromagnetic shower upon 
collision with the lead sheet. Electrons freed by the resulting ionisation of the gas 
are avalanched towards the anode wires where the resulting charge pulse capaci- 
tively induces a signal on nearby cathode pads (each measuring 30 mm x 30 mm). 
Thus the pads measure the position of the e.m. showers whilst the wires, which 
have a fast readout, can be used in the trigger. The 45 layers of a module are split 
into three 'stacks'. The first stack contains the first 10 layers and is 4 radiation 
lengths thick, the second contains the next 23 and is 9 radiation lengths thick, 
the third contains the final 12 and is also 9 radiation lengths thick since its lead 
sheets are double the thickness (4 mm as opposed to 2 mm). Cathode pads from 
consecutive layers are grouped and connected internally into 'towers', which are 
skewed such as to point to the IP. Energy deposits from all pads belonging to a 
particular stack within a tower are summed in readout. This provides three energy 
readings per tower, the magnitudes of which (absolute and relative) depend on the 
incoming particle ID and energy. 

In excess of 73,000 towers provide the high granularity required for good e — n 
separation in jets. Solid angle coverage is 3.97T with dead regions due to module 
boundaries representing 2% of the barrel surface and 6% of the endcap surface. 
Endcap and barrel modules are offset from each other by 15° to prevent coincidence 
of boundaries. Calibration over the full energy range is performed using Bhabha, 
tt and 77 events, and also independently of any events using radioactive sources 
mounted close to the gas inlets and injected periodically into the gas system. The 
fractional energy resolution, AE/E, for electrons and photons is (18 GeV 5 /y/~E)%. 
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Figure 3.10: Portion of a separated ECAL layer. 
The spatial resolution is (6.8 mm 



3.8 The Hadron CALorimeter (HCAL) 

The HCAL forms the final barrier to particles travelling through the detector, and 
is designed to stop hadrons and provide a measure of their energy. Muons typically 
penetrate its full thickness, and their departure from the detector is tagged by muon 
chambers surrounding the main body of the HCAL. A particle at normal incidence 
encounters 1.2 m of iron, equivalent to 7.2 interaction lengths for hadrons. The 
HCAL also serves as the main support structure of ALEPH, and as the return 
path for the magnetic field. 

It consists of a barrel of twelve modules and two endcaps of six modules each. A 
module consists of 22 consecutive layers of 5 cm-thick iron and streamer tubes, and 
is typically ~7m long depending on position. Each has a final layer of 10 cm-thick 
iron. The tubes consist of 8 cm-wide 'comb' profile PVC extrusions contained in 
plastic boxes. Nine 'teeth' form eight cells, each containing an anode wire running 
the length of the tube in a gas mixture of 22.5% Ar to 47.5% CO2 to 30% isobutane. 



A photograph is shown in Figure 3.11 . The internal surfaces of the cells are coated 
in graphite paint. Showers from collisions in the iron layers create ionisation in the 
tubes, and the resulting electrons avalanche towards the anode wires. Signals are 
induced on cathodes on both sides of the tube. The upper (open) side supports 
copper pads, the lower supports aluminium strips which run the length of the tube. 
Like the ECAL, the pads are grouped into towers pointing to the IP. There are 
just under 2,700 towers in all, and the angular (A0 x Ad) size they present to the 
IP ranges from 3.7° x 3.0° at 9 = 90°, to 15.0° x 2.5° at 9 = 6°. The aluminium 
strips are 4 mm wide. Their output is binary, indicating whether or not a tube has 
been fired at least once. They provide a detailed two-dimensional view of shower 
propagation within the detector which is important for muon identification. 

Surrounding the iron structure are two double layers of tubes, 50 cm apart 
around the barrel, 40 cm apart at the endcaps. These are the muon chambers. 
Here the pad layers are replaced by additional strip layers which run orthogonally 
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Figure 3.11: Photograph of a streamer tube. 

to the tubes, enabling three-dimensional hits to be obtained for charged particles 
leaving the detector. Together with information from the main bulk of the detector 
this provides powerful muon identification. 

The energy resolution of the HCAL for pions is AE/E = (84GeVVv / ^)%. 

3.9 The Luminosity CALorimeter (LCAL) 

The LCAL is a lead and wire chamber sampling calorimeter which provides the 
main luminosity measurement for ALEPH. An accurate measurement of the lu- 
minosity is essential in order to obtain reaction cross-sections from event-rates. 
The measurement is obtained from the event-rate of the QED t-channel Bhabha 
process, which is well understood. The differential cross-section is proportional 
to # -4 , and receives corrections due to interference from the s-channel process, 
e + e~— >Z/7— >e + e~, which decrease with 9. Thus in order to obtain the high event- 
rates needed for an accurate measurement of the luminosity, and to make insignif- 
icant the s-channel corrections^, the LCAL sits at very low 9. 

The LCAL consists of two annuli about the beam-pipe, 25 radiation lengths 
thick, beginning at z = ±263 cm. Each covers the range 44 mr to 160 mr and is 
formed by two semi-annular modules with 38 sampling layers each. The modules are 
almost identical to ECAL modules, with layers consisting of lead and proportional 
wire chambers (see Figure |3.10p . Like the ECAL, the layers are grouped into three 
stacks, and pads are grouped into towers which point to the IP. There are 384 
towers per module. 

Bhabha events are tagged by requiring back-to-back hits in both units above 
a threshold energy. The integrated luminosity is obtained by dividing the num- 
ber of events seen by the theoretical cross-section multiplied by the experimental 
efficiency An accuracy of the order of 0.5% is obtained. 



4 necessary since they depend on the properties of the Z which were not accurately known at 
the start of LEP. 
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3.10 The trigger system. 

Necessary to the operation of ALEPH is a system that will judge when an event 
of note (e + e~ annihilation, Bhabha or two-photon) has taken place and will sub- 
sequently initiate the readout of the subdetectors. This is the trigger system. It 
allows the TPC to exist in its insensitive state ('gate closed', see Section p.5| ) for 
much of the running time without loss of efficiency It also reduces the amount of 
uninteresting data written to tape, and the dead-time of the detector that results 
from readout. 

The trigger decision is split into three levels, each based on increasingly more 
complex information. The Level-1 trigger reaches a conclusion within 5 /xs, less 
than the bunch crossing time. It requires a good charged track in the ITC and/or 
information from the calorimeters indicative of a particle deposit. Should the 
decision be 'yes', the TPC gate is held open and the Level-2 decision is processed 
within 50 fis. This is essentially a recalculation of the Level-1 decision based now 
on three-dimensional information from the TPC. In the case of a 'yes-decision' 
readout of the whole detector is initiated, otherwise data-acquisition is halted and 
reset. 

During LEP1 the Level-1 trigger performed much better than was anticipated, 
making Levels-2 and 3 relatively unimportant as they vetoed only a small fraction 
of events. Since running at LEP2 energies however, noisier beam conditions have 
led to a large increase in the number of unwanted background 'events'. As a result 
the importance of Level-2 has risen greatly, now reducing the event-rate by in excess 
of 50%. Level-3 is an offline trigger, remains unimportant and is not described here. 

3.11 Event reconstruction 

Event reconstruction refers to the processing that is performed on the raw data 
output from the detector in order to correlate the information and produce final 
data that is more representative of what actually occurred in the event. Most 
important from the point of view of the analysis described in this thesis, is the track 
reconstruction. The tracks are the reconstructed helical trajectories of the charged 
particles and are formed by connecting the hits produced in the tracking chambers. 
The reconstruction algorithm starts in the TPC, linking nearby hits to form track 
segments. Segments which are compatible with a single helical trajectory with an 
axis parallel to the z-axis are connected to make tracks. At least four TPC hits are 
required for a track. These are then extrapolated down into the ITC and VDET 
where compatible hits are added. If there are ITC hits left unassigned to a track 
then the algorithm attempts in a similar way to form these into tracks (which will 
then not have assigned TPC hits). These tracks are required to have a minimum 
of four ITC hits. The tracks are parameterised by a set of five helix parameters. 
These are: the inverse radius of a track's circular projection in the xy plane, 1/R; 
the tangent of its 'dip angle' (equal to the ratio of its longitudinal and transverse 
momentum components), tan(A); its distance of closest approach to the z-axis, do; 
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and the ratio of its y and x momentum components and its z coordinate at that 
point, 0o and zo respectively. A final track fit is performed using Kalman filter 
techniques [^J. 

Once an event is reconstructed, the energy flow algorithm is run. The purpose 
of this is principally to create associations between charged tracks and calorimeter 
deposits, enabling the identification of particles and an improvement in the overall 
energy resolution, and is described in detail in |H| and p8 |. Unfortunately energy 



flow was not written with charged particle decay in the tracking volume in mind. 
It disregards tracks with a do greater than 2 cm or a |zo| greater than 10 cm, and 
so is of limited use in this analysis. 



Chapter 4 
The Search 



4.1 Introduction 

This thesis details a search for the supersymmetric process by which (lightest) neu- 
tralinos are pair-produced in e + e~ collisions, each then independently and promptly 
decaying to a slepton plus corresponding lepton, the sleptons then travelling a 
measurable distance in the detector before themselves decaying to a lepton plus 
gravitino. 

e + e"^n^ 11 11 -> Gil Gil 

The final state is thus four leptons (two from the IP, two from kinks or with 
large doQ) plus missing energy (since the gravitinos will not be detected). Each 
neutralino decays independently giving six distinct channels, ee, p,£i, ff , ep,, ef 
and p,f. An example event is shown in Figure |4.1| . This process can contribute 
to the exclusion of the GMSB parameter space (or potential for GMSB discovery) 
since, for a significant portion of the kinematically allowed space, the cross-section 
for the neutralino is greater than that for the slepton, and so when m x > mj, 
slepton production via neutralino decay can dominate over direct production. 

The ultimate goal of the search was to discover evidence for supersymmetry. 
In the absence of any such evidence, it was to formulate model-independent cross- 
section limits for the process and to exclude as much of the kinematically allowed 
parameter space as possible under a minimal GMSB model. 

The search was performed using a total of 628 pb _1 of LEP2 data taken over 
a range of centre-of-mass energies from 189 to 208 GeV, recorded between 1998 



and 2000 (see Table f4.1|) . The data was stored on 'DST's. These contain all 
the information about the reconstructed ALEPH events (unlike the usual format 
of MINIs), but events must fulfil some basic criteria to be written to DST from 
'POT' (essentially the same but containing all events). An event must have 

• more than 3 GeV of charged energy coming from within 5 cm in do and 20 cm 
in z 



x See Section 3.11 
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Energy (GeV) 


Year 


Integrated Luminosity (pb : ) 


188.6 


1998 


174.2 


191.6 


1999 


29.0 


195.5 


1999 


79.9 


199.5 


1999 


86.3 


201.6 


1999 


41.9 


199.8 


2000 


0.8 


201.8 


2000 


0.7 


202.7 


2000 


1.8 


203.8 


2000 


8.3 


205.0 


2000 


71.6 


206.3 


2000 


65.5 


206.6 


2000 


60.8 


208.0 


2000 


7.3 


Total = 




628.0 



Table 4.1: Integrated luminosities recorded by ALEPH at each centre-of-mass en- 
ergy achieved by LEP between 1998 and 2000. Year 2000 energies formed a contin- 
uum and have been placed in 1 GeV bins about integer values. The energy shown 
is the luminosity weighted average. 

• OR 2 or more charged tracks with do and zq within this range 

• OR total ECAL wire energy greater than 15 GeV and ECAL \t \ < 500 ns (t 
is the difference between the time of the event as recorded by the ECAL wires 
and the time of the nearest LEP bunch crossing: events due to conventional 
e + e~ collision typically have \t \ < 100ns) 

• OR 1 or more photon candidates and no charged tracks. 

A method was required to distinguish events from the process being searched 
for ('signal' events) from events due to known Standard Model processes ('back- 
ground' events). The technique used is called Monte Carlo analysis, and requires 
the simulation of both the signal and the background. The simulated events are 
called Monte Carlo data. The method for distinguishing signal from background 
can then be developed using the Monte Carlo data, and subsequently applied to 
the real data. The results can be compared to the results that would be expected in 
the absence of the signal as given by the background Monte Carlo, and the degree 
of evidence for the signal process can be quantified. 



4.2 Monte Carlo simulation 



4.2.1 Signal 



Signal Monte Carlo data was produced using SUSYGEN 2.2 [|3T| before subsequent 
processing by GALEPH p2| (the ALEPH detector simulation program based on 
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GEANT 0,0) and JULIA [|3|] (the ALEPH event reconstruction program used 
to reconstruct both Monte Carlo and real data events). The sleptons were set 
stable at generator level, then tracked and decayed by GALEPHQ. 

The characteristics of a signal event as observed in the detector depend only 
on the neutralino mass, and the slepton masses and lifetimes. GMSB dictates 
that either the sleptons are degenerate in mass, and therefore share the role of 
the NLSP (the co-NLSP scenario), or that the stau is lighter and so is the sole 
NLSP (the stau-NLSP scenario). In the first case the mass degeneracy also leads 
to lifetime degeneracy, and so all three sleptons are described by a single mass and 
a single lifetime. In the second case the mass difference between the stau and the 
selectron/smuon becomes an extra free parameter. But it is of little importance 
since if the stau is significantly lighter than the selectron/smuon then the neutralino 
will decay almost exclusively to it, and so no channel other than the ff is of 
importance and the mass of the selectron/smuon is irrelevant. Thus although 
GMSB parameter space is six-dimensional, the working parameter space for this 
search can be reduced to three dimensions, (m x , tj), where I can be e, fi or r in 
the case of degeneracy, and r otherwise. These are of course not the only relevant 
parameters. The neutralino production cross-section and branching ratios are of 
critical importance to the success of the search since they define the event rates 
for each channel. But for the purposes of signal Monte Carlo generation and the 
formulation of the selections they are not important. 

For the purposes of Monte Carlo generation it is necessary to identify the ap- 
proximate region of parameter space in which the selections are likely to be sensi- 
tive. In particular, slepton lifetimes should be chosen such that a good mapping 
of the variation of the selection efficiency with slepton lifetime is obtained over the 
region of interest. Thus it is natural to convert slepton lifetime to slepton decay 
length, since this is the principal feature of the signal events, and thus will have 
a far more direct relationship to the selection efficiencies. The natural working 
parameter space is then (m x ,mj, dj), where dj is the slepton decay length in the 
detector. At the time of Monte Carlo generation, sensible values can be chosen 
for the slepton decay length based on the detector size and geometry, and then 
translated into slepton lifetimes ensuring a good efficiency mapping. The slepton 
decay length is a function of m x , mj, y/s as well as t~ v In general this function is 
not simple because the ionisation caused by the slepton as it passes through the 
material of the detector will cause it to slow down, and so the decay length is also a 
'function' of the detector set-up. However, it is not necessary that the translation 



2 Some modification to G ALEPH was needed for the tracking of sleptons. This was performed 
by Clemens Mannert. 
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from decay length to lifetime be exact, and the following approximation is used: 

d i 



where, 



7/3 = Piab/mj 



Piab = J E 2 rest - mj + {E rest V x / m x) 



2 

i 



Erest = (m 2 + mf - mf) /2m 



X > 



PX = y E beam ~ m \ ■ 

Here 7 and (3 pertain to the slepton in the lab frame, Pi a b is the (initial) mo- 
mentum of the slepton in the lab frame if it is emitted at 90° to the lab direction 
in the neutralino rest-frame, E rest is the energy of the slepton in the neutralino 
rest-frame, p x is the momentum of the neutralino in the lab frame, E^ eam is the 
LEP beam energy and m\ is the mass of the relevant lepton. In the case of a mixed 
channel (i.e. the two sleptons have different flavours) the lifetime is calculated with 
both lepton masses and the average is taken. 

Initial work showed that for a given decay length, the selection efficiencies were 
only weak functions of the neutralino and slepton masses. Thus Monte Carlo data 
was generated to focus on the variation of selection efficiency with decay-length 
more than on the masses. Eight decay-lengths were chosen, and seven points in 
(m x , mj) space were generated for each of them, apart from the two highest decay 
lengths for which an extra two points in mass space were used (since some initially 
surprising results were obtained for long decay lengths which will be described in 
Section |7.2p , making a total of sixty points. This was done for each channel at 
both the lowest and the highest LEP centre-of-mass energies under consideration: 
189 GeV and 208 GeV. An exception was the ee channel at 189 GeV which was 
used as a test case, and so 21 mass points were generated per decay length. 500 
events were generated at each point. The slepton masses ranged from 67 GeV to 
102 GeV and the neutralino masses from 69 GeV to 103 GeV with never less than 
a 1 GeV mass difference between them. The decay lengths ranged from 2 mm to 
20 m. The exact points at which Monte Carlo data was generated for each channel 
at each energy are given in Appendix |A]. Slepton mass degeneracy was always 
assumed in the generation of the mixed channels (eft, ef, /if). 

All the signal Monte Carlo data was run through the POT — > DST selection 



described in Section 4.1, and events failing were appropriately tagged and not 



included in any signal efficiencies. 



4.2.2 Background 



Due to the nature of the signal a large array of Standard Model processes form 
significant backgrounds to the search. Only processes which are very easily cut 
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against, or which have a very low cross-section can be disregarded as sources of 
background. A list of the backgrounds analysed, together with their cross-sections, 
the number of events that were analysed, and the generators used for their produc- 
tion, is given in Table |4.2j . At 189 GeV (at which 174 pb -1 of data was collected) 
roughly 20 times the number expected in data were analysed for each process. 
At 208 GeV (at which only 7.3 pb -1 of data was collected) roughly 200 times the 
number expected in data were analysed for each process. These factors were not 
achieved for the two-photon to tau and up/down quark processes however, since 
limited Monte Carlo samples were available. 



Process 


At 189 GeV 


At 208 GeV 


Generator 






er(pb) 


No. analysed 


a(pb) 


No. analysed 






qq 


99.4 


340,000 


78.8 


120,000 


KORALZ 4.2 




TT 


8.21 


30,000 


6.61 


10,000 


KORALZ 4.2 


|d 


WW 


16.6 


60,000 


17.5 


26,000 


KORALW 1.21 


Ann 


zz 


2.76 


10,000 


2.79 


4,000 


PYTHIA 5.7 


Zee 


99.1 


350,000 


98.9 


144,000 


PYTHIA 5.7 


o 


77 — > t 


431 


800,000 


461 


600,000 


PHOT02 


+j 
o 


77 — > ud 


487 


400,000 


487 


250,000 


PHOT02 


1 


77 — > c 


93.2 


160,000 


94.4 


50,000 


PHOT02 


EH 


77 — > s 


23.9 


40,000 


23.7 


17,000 


PHOT02 



Table 4.2: The Standard Model background processes that were analysed as back- 
grounds to the search, a is the production cross-section. The final states of annhila- 
tion processes are created via the annihilation of the electron and positron, whereas 
the final states of two-photon processes are created via the collision of two photons 
which originated from the electron and positron. For information about KORALZ, 



KORALW, PYTHIA and PHOT02 see [36], |1, [38| and [39| respectively 



The generated events were processed by GALEPH and JULIA in the same way 
as for signa|]. Kinematic cuts were applied in the generation of the two-photon 
Monte Carlo data. For all the stated two-photon processes, events were required 
to have a final state invariant mass greater than 2.5 GeV. For the two-photon 
to tau process, final states were required to have at least 0.15 GeV of transverse 
momentum. For the two-photon to quark processes a scattering angle of at least 
5 mr was required for the electron or positron. 

The samples were split into two equal halves, one to act as a guide in formulating 
the selections, the other to give the final background expectation once the selections 



3 This was performed by other ALEPH members, generating Monte Carlo data for general 
ALEPH use. 
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were finalised (so as to minimise bias in the background expectation). Exceptions 
were the two-photon to quark processes where the high cross-section and limited 
Monte Carlo samples meant that the same sample had to be used for both. 
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Chapter 5 
Event Analysis 



5.1 Introduction 



The technique for the selection of candidate signal events is a two-step process. 
Firstly, using the ALPHA |gDJ package which facilitates access to the full event 
information, preliminary data candidates are identified using some simple preselec- 
tion. For those events satisfying the preselection, variables are calculated that are 
designed to be sensitive to the differences between signal and background. These 
are written to HBOOK ntuples (essentially tables of data in which rows corre- 
spond to successive events and columns correspond to variables) forming a greatly 
reduced data set. The second stage is the application of a selection, consisting of 
a set of cuts on the variable values, to determine the final signal candidates. This 
is the technique used by most searches for physics beyond the Standard Model at 
LEP. 

Many, if not all, of the variables should pertain to the event as a whole, such 
as the number of charged tracks, total energy and invariant mass. In order to 
calculate such 'global' variables the relationships between the various reconstructed 
detector objects such as tracks and calorimeter deposits should be known. This 
is the purpose of the energy flow algorithm. But energy flow has only limited 
abilities to make sense of processes that are happening away from the interaction 
point. It ignores tracks with a do greater than 2 cm or a z greater than 10 cm, and 
since JULIA does not reconstruct the slepton decay vertex, it does not recognise 
the relationship between the reconstructed slepton track (if there is one) and the 
resulting lepton track. Thus energy flow will only yield correct results if the slepton 
decay length is short enough for the resulting lepton to satisfy the cuts on d and 
z , or long enough for the slepton to reach the calorimeters before decay; and so it 
is of limited use in this analysis. 

Even if energy flow did treat charged particle decay in the tracking volume 
correctly, the signal would still not be well-defined by global variables in the case 
of the ff channel (and the ef and jlf channels to a lesser extent) because of 
the larger number of invisible particles (up to eight neutrinos on top of the two 
gravitinos) and the larger range in the charged track multiplicity (from multi-prong 
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tau decays). The selection then relies primarily on the presence of high-do tracks: 
that being the primary feature of the signal. Several other processes can generate 
high-do tracks though, such as nuclear interactions, ECAL splash-backs, photon 
conversions, multiple scattering and cosmic rays. The frequency with which these 
processes occur means that their rejection efficiency must be close to 100% if they 
are not to swamp any signal. 

After the event preselection is described in Section |5.2| , this chapter goes on to 
detail the procedures that are used to discriminate between high-do tracks due to 
the signal process, and those due to background. Sections fy| and [5]4] describe the 
identification of ECAL splash-backs and nuclear interactions respectively. Then 
Section ^75] describes the procedure that is used to confirm (or reject) that a high- 
do track is compatible with the hypothesis of originating from slepton decay. The 
final cut-based selections are described in Chapter ||. In order to facilitate the 
descriptions of the procedures described in this chapter a terminology is used which 



is explained below and in Figure |5.1| . The figure also introduces some points of 
interest concerning tracks reconstructed from particles that have not been produced 
at the IP. There are three terms which are important to understand: 

The post(pre)-do trajectory of a particle. This is the section of the par- 
ticle's trajectory after it has passed (before it has reached) its point of closest 
approach to the beam-axis. 

A post(pre)-d track. This is a track that has been formed by hits created 
on the post(pre)-do portion of a particle's trajectory. Both this and the term above 
refer to truths about the event that are not known a priori from the reconstructed 
information. 

The post(pre)-d section of a track. This is the section of a track which is 
post(pre)-do according to the track direction (assigned by JULIA, and which may 
be incorrect). The track will only have assigned hits on its post-do section. If it is a 
pre-do track then its pre-do section will probably follow closely the post-do section 
of a post-do track (see Figure |5Tl| ). 

Although this terminology might seem confusing, it is an unfortunate reflection 
of the subtleties involved in the relationship between the reconstructed and 'truth' 
information of tracks not originating from the IP. 

Many of the procedures involve the fitting of tracks to vertices, either to test 
whether a group of tracks form a new vertex, or whether they fit to an existing 
vertex. This is performed using the ALEPH YTOP package described in and 
43|j . It provides a y 2 for each fit performed, which is normalised to the number of 



degrees of freedom. 



5.2 Preselection 

In order to speed up analysis, events were first subjected to some very basic pre- 
selection aimed at fast rejection of events that were clearly not signal candidates. 
The event was rejected if 

• The scalar sum of the momenta of all tracks was less than 4 x [Elep /188.6) GeV 
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Figure 5.1: The reconstructed track(s) of a particle that has been produced in the 
TPC, moving towards the beam-axis, through the ITC, then back out into the TPC 
and on to the ECAL. The do-point is the point at which the track is closest to the 
beam-axis. Before it reaches the do-point it is travelling along its pre-do trajectory, 
afterwards, along its post-do trajectory. Since the idea that particles are produced 
at the IP is a paradigm within JULIA, it assumes that the particles responsible for 
all tracks were moving away from the beam-axis and did not originate significantly 
before their do-point. For this reason it will not form tracks from hits that would 
lie on both sides of the do-point; and so the particle trajectory shown here is 
reconstructed as two tracks: one pre-do track reconstructed from the pre-do hits, 
and one post-d track reconstructed from the post-d hits. JULIA incorrectly 
assigns the directions of pre-d tracks as being radially outwards, and so both their 
momentum 3-vector and charge have the wrong sign. Tracks that have no TPC 
hits will be referred to as ITC tracks, tracks that do have TPC hits will be referred 
to as TPC tracks and may or may not have ITC hits. Either type of track may or 
may not have VDET hits. 
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(to reduce the number of two-photon events) 

• The total number of tracks was less than 3 or greater than 30 (to reduce both 
two-photon and qq events) 

• The total number of tracks was less than 4 and the two highest momentum 
tracks had momenta greater than half the beam energy and the angle between 
them was greater than 179° (to reduce e + e~ and events) 

• The total number of energy flow charged tracks was greater than 11 (to 
further reduce qq events) 

• The energy flow event energy was less than 7 x (Ex^p/188.6) GeV and the 
energy flow event momentum was less than 1 x (.Elsp/188.6) GeV (to further 
reduce two-photon events) 



5.3 Rejection of ECAL splash-backs 

The electromagnetic shower created when a charged particle enters the electro- 
magnetic calorimeter typically propagates in the direction of the incoming particle. 
Some of the electrons and positrons in the shower can however, scatter back into 
the TPC where they can be reconstructed as tracks. This is known as an 'ECAL 
splash-back'. 




Figure 5.2: An example of an ECAL splash-back in a Monte Carlo rr event at 
189 GeV. Two electrons have escaped the ECAL to be reconstructed as high do 
tracks in the TPC, one arcs through the ITC and passes out the other side to form 
a third track. 
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They are characterised by a high-do track 'entering' the ECAL (in fact the 
particle responsible is leaving but JULIA assumes all tracks move radially away 
from the beam-axis: the track is pre-do) close to another track that has higher 
momentum, with both forming a common vertex in the region of the ECAL (see 
Figure 1572] for an example). The only circumstance under which this is likely to 



occur in a signal event is when the slepton decays inside the ECAL and the resulting 
lepton is ejected back into the TPC. Since this analysis is only intended to be 
sensitive to slepton decay in the tracking volume then the possibility of rejecting 
signal events by confusing slepton decays with ECAL splash-backs is not regarded 
a problem. 

The track that results from the splash-back will be referred to as the output 
track, the track causing the splash-back will be referred to as the input track. The 
conditions that must be satisfied by both a candidate output and input track for 
the output track to be tagged as a splash-back are as follows: 

• The output track must have TPC hits and momentum, p ou t pu t < 
5 x (E LEP /188.Q) GeV and d > 0.5 cm or z > 2 cm. 

• The input track must also have TPC hits and have Pi npu t > 0.5 GeV and 

Pinput Poutput- 

• The points at which the output and input tracks leave the TPC volume must 
be within 50 cm of each other. 

• The tracks must form a common vertex in the region 

165 < r < 230 cm and \z\ < 300 cm, 
or 200 < \z\ < 300 cm and r < 230 cm 

(these roughly correspond to the vicinity of the ECAL barrel and endcaps 
respectively). 



5.4 Rejection of nuclear interactions 

Hadrons interacting with nuclei in the structure of the detector can cause hadronic 
showers leading to multiple high-do tracks. These are characterised by many tracks 



originating from a common vertex which is displaced from the IP (see Figure [5T3 
for an example). Like splash-backs these have distinct characteristics that can be 
tagged. Care must be taken though to avoid confusing a f — > r mu i ti _ prong decay 
vertex with a nuclear interaction. 

Clearly an algorithm that can find multiple tracks originating from a common 
vertex is required. Initially this was done by testing all track pairs for compatibility 
with a common vertex, and then attempting to merge these 2-track vertices. This 
was very CPU intensive however, and a method which used the fitting routine less 
often was sought. The procedure settled upon relied on starting with a single high- 
do track, and trying to add further tracks to it to build up a multi-track vertex. 
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Figure 5.3: An example of a nuclear interaction in a Monte Carlo rr event at 
189 GeV. A pion from the decay of a tau has interacted with a nucleus in the inner 
TPC wall, giving a shower of high-do tracks, mainly protons. 



Once all the tracks that could be added had been (which may be none at all) , and 
if there were one or more other high-do tracks in the event which had not been 
included in the vertex, then one of these was used in the same way as the first 
track to try to build a new vertex. Unlike the first method tried, this method 
required making a choice as to the order in which the tracks should be selected so 
as to give the best chance of all tracks from a given vertex actually being included. 
Using tracks in order of increasing fractional momentum error was found to be 



successful. The exact procedure used is shown in the flow diagram of Figure |x4 



In this diagram 'most reliable' means smallest fractional momentum error. A track 
passes 'basic cuts' if it has not been tagged as an ECAL splash-back, has TPC hits, 
has a momentum smaller than half the beam energy and greater than 0.1 GeV, and 
a do> 0.15 cm or a zq> 3 cm. A track is 'hit-compatible' with a vertex if at least 
all bar one of its hits lie on one side of the vertex, so that it is compatible with the 
hypothesis that it was produced at that vertex. 

Once one multi-track vertex has been found with this procedure, then if each 
new vertex after the first is less than 3 cm from, or shares more than 50% of its 
tracks with, a previously found vertex then an attempt is made to merge the two. 
If this results in a new vertex with a \ 2 (given by YTOP) less than 70 then the 
merge is considered successful, and the two vertices are replaced with the single 
merged vertex. 

Once the vertices have been finalised, certain properties are calculated for each, 
such as the summed momentum and charge of the constituent tracks. Care was 
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Track 1 = most reliable track 



NO 



Is there a less reliable 
track than track 1 ? 



YES 



YES 



NO 



Exit 



1_ 



Is track 1 already included in a vertex? 



NO 



Track 1 = next 
most reliable track 
after existing track 1 



NO 


Does track 1 


YES 


< 


pass basic cuts? 


• 



Track 2 = next 
most reliable track 
after track 1 



Is there a less reliable 
track than track 2? 



YES 



Track 2 = next 
most reliable track 
after existing track 2 



Note: track 1 may 
now be included 
in a vertex 



YES 



NO 



YES 



Does track 2 pass basic cuts? 



Is track 1 now included in a vertex? 



Is track 2 compatible with 
current vertex coordinates? 



YES 



NO 



Does track 2 complete less 
than 3/10ths of a full spiral 
between its dO and the vertex? 



YES 



NO 



NO 



If both tracks have dO<lcm, 
is angle between them >25 degrees 
and the radius at which 
they cross >1 cm? 



YES 



Does track 2 fit to the 
current vertex giving a vertex 
chi-squared less than 70? 



YES 



NO 



NO 



Do they form a common 
vertex with a chi-squared 
less than 70? 



YES 



YES 



Does the inclusion of track 2 
cause the vertex to move such 
that one or more other tracks 
are no longer hit-compatible? 



NO 



Are both tracks hit- 
compatible with the vertex? 



NO • 




Track 2 is included in the vertex 



YES 



NO 



Tracks 1 and 2 form a new vertex 



Figure 5.4: Flow diagram showing the procedure used to find and tag nuclear 
interaction vertices. 
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taken however, to account for a potentially problematic feature of the reconstruc- 
tion algorithm. Since JULIA always assumes that particles move radially away 
from the beam axis, a particle that is produced far enough from the beam axis 
and moving back towards it can result in two reconstructed tracks (a pre-do and 
a post-do), as is shown in Figure |5.1|. In this case the track corresponding to the 



portion of the particle trajectory before it reaches its point of closest approach 
to the beam axis is assigned the wrong direction. Thus its momentum 3-vector 
and charge have the wrong sign, and this will corrupt the overall properties of any 
multi-track vertex in which it is included. Also, since the total number of tracks 
in the vertex is used to classify it, then a single particle giving two tracks would 
corrupt this as well. To correct this, any track whose TPC hits lie at a smaller 
radius than the vertex has its charge and momentum 3-vector reversed. Then the 
do, z and directions of these tracks are compared with any tracks which are part 
of the same vertex and for which the vertex lies on the pre-do portion. If the tracks 
differ in d and z by less than 30% of their respective average, and the angle be- 
tween them at their do-points is greater than 170°, then both are assumed to have 
originated from the same particle trajectory and the one with fewer TPC hits is 



removed from the vertex. As an example, if both the tracks shown in Figure pTL 



were found to originate from a nuclear interaction, then the pre-do track would be 
removed. 

Each vertex is then placed into one of 3 categories based on its properties. A 
vertex is labelled as a f — > T multi _ prong vertex if it satisfies all of the following 
conditions - 



• its charge is ±1 

• its largest two-track angle at the vertex is less than 20° 

• its momentum is > 4 x (E LEP /188.6) GeV and < 55 x (E LEP /188.6) GeV 

• its kink- angle is between 5° and 170° 

• its kink-angle is > 15° or its | cos(#)| is less than 0.6 

• its number of tracks is less than 6 



• its distance from the IP is greater than 0.3 cm 

• and its distance from the IP is greater than 1 cm or its \ 2 is l ess than 3. 

The charge of the vertex is the sum of the charges of all the constituent tracks. 
The momentum is the vector sum of the momentum 3- vectors of all the constituent 
tracks. The kink-angle is the angle between the vertex momentum 3-vector and 
the straight line from the IP to the vertex. The | cos(#)| of the vertex refers to the 

Vertex position, i.e. \z ver t ex / "\/ r vertex z vertex I ■ 

Although this is a reasonably stringent set of cuts, background events can 
still have vertices labelled as f — > r mu i ti _ prong since a charged particle emitting a 
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bremsstrahlung photon which quickly undergoes pair-conversion can create some- 
thing very similar. 

If a vertex fails this selection then it is classified as a nuclear interaction vertex 

if- 

• its radius is greater than 4.8 cm (a nuclear-interaction is highly unlikely inside 
the beam-pipe, so the minimum radius is set at the beam-pipe radius of 
5.3 cm, minus 5 mm leeway) 

• its largest two-track angle at the vertex is greater than 20° 

• and its number of tracks is greater than 2. 

If it fails this also then it is classified as 'unknown', and the constituent tracks are 
still considered valid. 

5.5 Confirmation of slepton-decay hypothesis 

Although this analysis attempts to reject high-do tracks from known sources such 
as those discussed in Sections [|]3] and |5]4], and these attempts are largely successful, 
they are not sufficient. The number of residual high-do tracks is still too many. 
Hence a procedure was developed, not to test whether a high-do track comes from 
a known background process, but to test whether it is compatible with coming 
from the signal process, i.e. the decay of a heavy charged particle. The hits of a 
candidate track are examined to identify the region of its length over which the 
responsible particle could have originated, and this is then used to calculate what 
can be expected as the detector response to the parent slepton. In the absence of 
such a response the track can be rejected. For instance, if a track's hits indicate that 
the responsible particle must have been produced far out in the TPC, the parent 
slepton should have its own reconstructed track which forms a good kink-vertex 
with the high-do track. In the absence of such a 2-track vertex the high-do track 
can be rejected. This section describes the exact procedure in detail. It should be 
noted that only tracks that have TPC hits, a momentum greater than 1 GeV, a do 
greater than 0.1 mm and a Xbs X 2 °f the track- fit to the beam-spot per degree 
of freedom) greater than 5 are considered as candidate high-d tracks. 

5.5.1 A particle's point of origin 

While a reconstructed track is a reliable indicator of a particle's trajectory, it does 
not tell us exactly where that particle originated or where it stopped/decayed. 
Nominally, the point of origin is considered to be the do-point. This is fine for 
tracks originating from the IP, since the distance scale associated with the uncer- 
tainty on the particle's exact point of origin will be at most the size of the beam 
spot, which will be several orders of magnitude lower than the track's radius of 
curvature. Thus the track is effectively straight over the relevant range, and the 
uncertainty on its point of origin does not translate into significant uncertainty on 
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its initial momentum 3-vector. For a track not originating from the IP no such 
simple assumption can be made, and only analysis of its hits can yield information 
as to its point of origin. This will take the form of a range of the track length over 
which the particle could have originated. 

The following procedure for finding this range assumes that the track in question 
is post-do (i.e. reconstructed from hits created as the particle moves away from 
the beam axis), and so the track direction and charge are correctly assigned by 
JULIA. Although the following procedure is not sensible for a pre-do track, it will 
be explained later that such a track is in fact dealt with in the correct way. 

The range over which the particle could have been produced is between the 
lowest-radius point at which there is evidence that the particle was following the 
track trajectory, and the highest-radius point at which there is evidence that it 
was not. To find evidence that the particle was not following the track trajectory, 
a scan is performed backwards along the track length, starting from the point at 
which it achieves its highest radius in the TPC, looking for missing hits. Nominally, 
a hit is expected at every TPC pad radius that the track crosses, every ITC wire 
radius that the track crosses while it is inside the ITC active volume, and every 
intersection that the track makes with a VDET wafer. There are several legitimate 
reasons why a track can be missing a hit however. The TPC sector boundaries 
form dead regions of the order of 1 cm wide, and in the ITC certain wires are not 
active. Also, even if the particle has produced a hit it may not be assigned to the 
respective track since the reconstruction algorithms can simply fail to assign the 
hits correctly, or may only have one hit to assign between two tracks (if they are 
closer than the minimum two-track resolution). Thus if the track fulfils any of the 
conditions given in Table 5T at a point where a hit is expected but not observed, 
the lack of a hit is not considered abnormal. 

If a point is encountered at which there is no hit and at which none of the 
conditions of Table |5J] are satisfied, then it is assumed that the responsible particle 
was not following the track trajectory at this point. Then if there are no hits at 
lower radii, the range for the particle point of origin is set between this point and 
the lowest-radius assigned hit. If there are one or more hits at lower radii, then the 
track information appears to be inconsistent with a genuine particle trajectory. In 
this case, if the point under consideration is in the VDET the track is rejected as 
bad. If it is in the ITC then the track is spared for now and the scan continues, 
but if another such point is encountered the track is rejected. If the point is in 
the TPC, then the location of the bad point is merely recorded. These bad TPC 
points will be referred to again later. 

If the scan reaches the track's do-point without the track being rejected or the 
particle range-of-origin being determined then the scan continues along the pre-do 
section of the track, as far as the highest radius that the pre-do section of the 
track reaches in the TPC if necessary. On the pre-do section the track will not 
have any assigned hits and so Table 5A_ governs completely whether the track is 
considered good at each point. If the scan reaches the end of the pre-d section 
without conclusion, then the track is rejected (since it is probably a cosmic muon 
that has passed right through the detector). 
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TPC 


• The track is < 1.5 cm in rip from a sector boundary 

• The track is < 1.5 cm in r<p and < 4 cm (10 cm) in z from another 
TPC (ITC) track on the same side (on either side) of the central 
membrane 

• The track is < 2 cm in rcf) and < 4 cm in z from an unassigned hit 

• The radius of the point in question is < 3 mm greater than the 
track's do 


ITC 


• The nearest wire has a hit 

• The nearest wire is dead 

• The track is < 1 mm in r<p from the cell edge and the next-nearest 
wire has a hit 

• The point in question is at a radius < 0.1mm above the track's do 


VDET 


• There is a hit within 0.5° in <fi 

• The track is within one sigma of the wafer edge 



Table 5.1: Conditions under which the lack of a hit is not considered abnormal in 
each tracking subdetector. 



5.5.2 Where is the slepton? 

Once a range for the particle's point of origin has been calculated, it is possible 
to translate this into a range of possible trajectories for the parent slepton, and to 
make the requirement that the event information supports one of these trajectories 
and thus the hypothesis that the candidate track originated from slepton decay. 
This section describes the procedure that is applied to do this. 

There are four distinct scenarios for the reconstructed detector response to 
the slepton depending on how far the slepton travels before decay and in what 
direction. If the slepton achieves a high enough radius before decay it can produce 
a reconstructed TPC track, which should form a well-reconstructed vertex with the 
lepton track. In the following discussion this will be referred to as a TPC-vertex. 
Failing this, it can produce an ITC track which will form a vertex with the lepton 
track that is not so well reconstructed in z due to the low z-resolution of the ITC. 
This will be referred to as an ITC-vertex. Failing this, if the slepton at least passes 
through some active tracking components before decay, then it should leave some 
hits which point to the lepton track. These will be referred to as pointing-hits. 
Finally there is the case that the slepton decays before creating any hits, and so 
the only sign of its existence is the do of the lepton track. Examples of each case are 



shown in Figure |5.5| . The range for the point of origin of the particle responsible 
for the candidate high-do track determines which of these scenarios are possible, 
and which are not. 

The relationship between the point of origin range and the possible scenarios is 
as follows. If, for at least part of the range, 
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c). d). 

Figure 5.5: Examples of the typical reconstructed detector response to the slepton 
and lepton for decreasing slepton decay lengths. Shown are (a) a TPC-vertex, (b) 
an ITC-vertex, (c) pointing-hits and (d) a simple high-do track from signal Monte 
Carlo data. All tracks and hits other than those of the slepton-lepton pair shown 
have been removed for clarity. 
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• the track radius (i.e. the distance of the track from the z-axis) is above that 
of the 4th TPC pad row, then a TPC-vertex is possible. 

• the track radius is above that of the 4th ITC wire layer and the track is either 
inside the ITC or the straight line from the track to the IP enters the ITC at 
a radius greater than this (and thus the slepton, under the assumption that 
its trajectory was straight, crossed at least 4 ITC wires), then an ITC- vertex 
is possible. 

• the straight line from the track to the IP passes through some tracking com- 
ponents, then pointing-hits are possible. 

• the straight line from the track to the IP passes through no tracking compo- 
nents, then no slepton response at all is possible. 

These possibilities are not mutually exclusive, they could all be open. However, 
of those that are open, it may be possible to exclude some based on likelihood. 
For instance, if the track radius is significantly above that of the 4th TPC pad 
row for the entire range, then anything other than a TPC-vertex is unlikely, and 
a TPC-vertex can be demanded. Thus if any of these conditions, with 3 cm added 
to the radius necessary for a TPC-vertex and 3 mm added to the radius necessary 
for an ITC-vertex, are true over the entire range, then the first of the respective 
possibilities becomes a demand. If the demand is not satisfied, the track is rejected. 
If it is found that no single response can be demanded in this way, then the track 
is tested for the possible slepton responses in the order TPC-vertex, ITC-vertex, 
pointing-hits. If the track does not satisfy any of the possible responses (bearing in 
mind that if the case that the slepton did not produce any hits has been identified 
as a possibility, then no response is possible), then it is rejected. 

The procedure involved in looking for TPC-vertices, ITC- vertices and pointing- 
hits will now be described. 

TPC-vertices 

The first step of testing a track for a TPC-vertex is to identify candidates for the 
slepton track. These are required to have TPC hits (by definition), a do less than 
2 cm and a zo less than 2.5 cm, and must not have been identified as tracks from 
a nuclear-interaction or an ECAL splash-back. Its dE/dx (if dE/dx information 
is available for the track) must be no less than 10 standard deviations below that 
expected for a singly-charged particle of mass 60 GeV (the approximate GMSB 
slepton mass limit), and the track should show evidence that the particle respon- 
sible ceased to follow the track trajectory before leaving the TPC (i.e. be missing 
at least one hit from its end without observable reason - see Table |5.1|) . If there 
are no candidate slepton-tracks, there can be no TPC-vertices. 

The high-do track is tested for a good TPC-vertex with each candidate slepton 
track. They must form a common vertex which lies on the post-do section of the 
slepton track and has a \ 2 (from YTOP) less than 30. The z-coordinates of the 
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Figure 5.6: An example of a TPC-vertex with a decay angle greater than 90°. 
The lepton has formed a pre-do and a post-do track. The arrows show the track 
directions as assigned by JULIA. 

two tracks at the vertex must differ by less than 7 cm, the angle between them at 
the vertex must be greater than 10°, and the last assigned hit of the slepton track 
must not be at a radius more than 1 cm above that of the vertex. 

The checks made on the high-do track's hits are more complex, since there 
are three possibilities for their locations. If the decay-angle (defined here as the 
angle between the slepton and lepton directions at the point of decay in the r0 
plane) is less than 90° then all the high-do track's hits lie at radii greater than the 
vertex radius, and both the track directions are correct and so the decay-angle will 
be calculated correctly (see Figure |^.a). If the decay-angle is greater than 90° 
and the high-do track is post-do, then again the track directions are correct and 
the decay-angle will be calculated correctly, but now the hits can lie both above 
and below the vertex radius so no constraints can be made on their positions (see 
Figure |5.6|) . However, if the decay-angle is greater than 90° and the high-do track 
is pre-do ( so & U its hits lie at a radius smaller than that of the vertex), then the 
decay-angle will be calculated incorrectly to be less than 90°. So if the calculated 
decay-angle is less than 90°, all the hits are required to lie either at radii above 
that of the vertex minus 1 cm, or at radii below that of the vertex plus 1 cm. If 
the calculated decay angle is greater than 90° then no constraints are made on the 
positions of the hits. Thus if a high-do track is found to have a TPC-vertex it 
can be tagged as post-do or pre-do based on the positions of its hits relative to the 
vertex. 

At this stage, the number of bad TPC-points for the track (see Section |5.5.1|) 
can be fixed. For pre-do tracks it is set as the number lying at least 6.5 cm (normal 
pad-row spacing plus 1 mm) below the radius of the vertex. For post-do tracks it is 
the number lying above the radius. The number of bad TPC points forms a track 
variable which can be used as an indicator of the reliability of the track, and is 
referred to in the selections described in the next chapter. 

ITC-vertices 

As with TPC-vertices, the first step towards finding ITC-vertices is to identify 
candidate slepton tracks. These are required to have no TPC hits (by definition), 
a do less than 1 cm (no constraint on zq due to the low ITC ^-resolution) and a 
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momentum greater than 1 GeV. 

The high-do track is then tested for a good ITC-vertex with each candidate 
slepton track. Because of the low ^-resolution of the ITC a full vertex fit is not 
performed. It is merely required either that the tracks' circular r<ft projections cross 
in the r<j) plane, in which case the crossing point closest to the beam axis is taken 
as the two-dimensional vertex, or that their distance of closest approach in r<ft is 
less than 5 mm, in which case the vertex position is the point of closest approach. 

The vertex radius is required to be no more than 5 mm below that of the 4th ITC 
wire layer. The upper limit on the radius is determined by a scan along the ITC 
track outwards into the TPC, looking for evidence that the particle responsible was 
not following the track trajectory. At successive pad-row radii a set of conditions 



similar to those summarised in the TPC section of Table 5.1 are evaluated, with 



the z requirements relaxedP]. The upper limit on the radius is then set at the radius 
of the second point at which the track is missing a hit without observable reason 
plus 5 mm, unless the first and second points are consecutive in which case the 
first is used. The first point is in general not used since the track errors in the 
TPC are much greater in the absence of assigned TPC hits, and so some leeway is 
appropriate. 

The ITC track is required not to have any assigned hits on wire-layers above 
the layer that is itself just above the vertex. The constraints on the TPC hits of 
the high-do track are the same as for TPC-vertices, and so if it passes it is known 
whether it is post-d or pre-d , and its number of bad TPC points is fixed in the 
same way. If it is found to be post-d with a decay-angle less than 90° then it is 
further required to have no ITC hits below the wire layer that itself is just below 
the vertex, and no VDET hits at all. 

While TPC-vertices are very rare in background events, ITC- vertices are less so 
because of the weaker constraints in the z direction, and because multiple scattering 
in the ITC outer wall and TPC inner wall can cause the ITC and TPC hits of a 
particle to be reconstructed as separate tracks, giving vertices very similar to those 
in signal. Some rejection of these is attained by the requirement that the decay- 
angle (corrected in the case of pre-d tracks) is greater than 15°. 



Pointing-hits 

The object here is to find a line of hits compatible with marking the slepton's 
trajectory from the IP to the high-do track. Although the slepton's trajectory is 
helical, it should have sufficient momentum that it can be considered straight over 
the distances involved. This is done to keep matters simple, but a small curvature 
of the slepton's path can be accommodated in the leeway used in the following 
procedure. 

The two endpoints of the particle's point-of-origin range are used to define the 



lr rhe separation in z for proximity to another track to be considered a valid reason for missing 
a hit is set at 10 cm for both ITC and TPC tracks, and no z information at all is used in judging 
the proximity of unassigned hits. 
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0-range in which the line of hits should lie. If the track lies at sufficiently high 
| cos(6 l )| over the whole range, then the straight line from the IP to the track at 
any point in the range will not cross any tracking components and it is therefore 
possible that the slepton did not produce any hits. The track is accepted without 
further requirements in this case. If this is not the case, the procedure continues 
to look for pointing-hits using hypothesised slepton trajectories. 

Firstly it is determined if a slepton hit is required in the first VDET layer. This 
is the case if the track radius is above the minimum VDET layer- 1 radius (VDET 
layers are not cylindrical and so do not sit at a single radius), and the track | cos(0)| 
is below the maximum VDET layer- 1 | cos(#)| for the entire range. In this case, 
unassigned layer-1 VDET hits are sought in the relevant (f> range. If none are found 
the track is rejected. If one or more are found, they dictate the possible trajectories 
for the parent slepton. These are straight lines from the IP with 0's equal to those 
of the VDET hits, and #'s such that they intersect the track. If they also intersect 
the second VDET layer, an unassigned hit is required here also, with a within 1.2° 
of that of the first VDET hit (i.e. no more than ~ 2 mm from the slepton trajectory 
given by the first hit). Trajectories for which a second VDET hit is expected but 
not observed are rejected. If a slepton hit in the VDET is only required for part of 
the track's point-of-origin range, then VDET hits in the relevant range are still 
sought, but the track is not rejected if none are found. For the part of the range 
where slepton hits are not required (which may be all of the range), hypothesis 
slepton trajectories are taken from the <f>'s of the ITC wires in the first layer that 
lie in the relevant <ft range. For each trajectory which crosses at least one ITC 
wire-layer, ITC information is required to support the hypothesis that the slepton 
followed this trajectory. This is tested by applying the conditions summarised in 
the ITC section of Table |5.1| at each wire-layer radius. If a slepton trajectory is 
found not to be valid in this way, it is rejected. If the remaining number of slepton 
trajectories is non-zero, the track is considered to have pointing-hits. 



A note on pre-do tracks 



It was mentioned in Section [5.5.1 that the procedure for calculating the point- 



of-origin range assumes the track is post-d , and is nonsensical for pre-d tracks. 
However, the ordering in which the procedures are applied in the computer code is 
not exactly as is implied in this chapter. In particular, the procedure for finding 
TPC-vertices precludes the possibility that a TPC-vertex could be found which is 
incompatible with the relevant track's point of origin range. So the first thing that 
is actually done is that a TPC-vertex is sought. If one is found then the track is 
accepted without the particle point-of-origin range ever being calculated. While 
this makes sense from a pragmatic viewpoint since it speeds up the analysis, it also 
solves the problem of pre-do tracks. Since a pre-do TPC track can only be formed 
if the particle point-of-origin is at a radius greater than that of the 4th TPC pad 
row (the particle will need to create at least 4 TPC hits before the do-point), it is 
reasonably safe to say that the slepton must also have crossed 4 TPC pad rows, 
and so also have a reconstructed TPC track, and so the track should have a TPC- 
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vertex. The unconditional acceptance of tracks with TPC-vertices therefore means 
that pre-do tracks are treated correctly. 

5.5.3 The track variable, 

All tracks which are tested for compatibility with the hypothesis that they origi- 
nated from slepton decay using the procedures described in this chapter, are as- 
signed an integer value to indicate the level and type of success. This value is the 
track variable, Jij. 

An ftj of indicates that the track is not compatible with the hypothesis, and 
such tracks are rejected as good high-do tracks. A track has an Oij of 5 if it has a 
TPC-vertex, an Jlj of 4 if it has an ITC-vertex, an %j of 3 if it has pointing-hits, 
and an Oij of 1 or 2 if the track is compatible with the hypothesis that the slepton 
did not produce any hits. A value of 2 is given in the case that the track has 
sufficiently small do that it is possible that the slepton decayed before reaching a 
high enough radius to produce a hit. A value of 1 is assigned when the slepton did 
reach a high enough radius, but had a | cos(0) | high enough that it still managed not 
to pass through any tracking components before decay. This distinction is drawn 
since an Oij of 1 is the most unlikely of the values under the signal hypothesis, but 
this relative probability is not mirrored in the background, where detector effects, 
cosmic muons and beam-gas events mean that an "Jif of 1 is not so improbable. 

Examples of tracks with %fs from 2 to 5 are shown in Figure |5.5| , and an 
example of a track with an of 1 is shown in Figure [577. 
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Figure 5.7: An example of a track with Jij = 1 from signal Monte Carlo, viewed 
in (a) the r<p plane and (b) the pz plane (p is equal in magnitude to r, but with 
a sign that makes it negative for a certain half of the range which is chosen to 
make the event as clear as possible). Although the do of the track shows that the 
slepton reached a radius above at least one VDET layer before decaying, it did so 
at such a high | cos(#) | that it did not produce hits. 



Chapter 6 
Selections 



This chapter describes the sets of cuts that are made on the variables calculated 
from the full event information that form the final discrimination between signal 
and background. There are two sets of six selections. The six selections of a given 
set are aimed at selecting with high efficiency events corresponding to each of the six 
possible final states (ee, \x\x, ff, ejl, ef and /if), while minimising the acceptance for 
background events. The difference between the two sets is only that the three most 
critical cuts are tuned under different conditions, such that one set gives optimal 
sensitivity in the case of low slepton decay-length, and the other in the case of 
high slepton decay-length. As such the two sets of six selections will be called 
the low and high decay-length optimisations. Both sets are applied independently 
to the real data and the Monte Carlo data, and so there are two sets of results 
described in Chapter [7[ More details on the "hows and whys" of optimisation will 
be given in Section |6.4| . The exact differences between the optimisations are given 
in Section |6T| 

Each selection is split into track and global cuts. Track cuts are cuts applied 
to variables that pertain to a given track within an event and are used to select 
good high-do tracks. Global cuts are cuts on variables that pertain to the event as 
a whole, such as the total energy or the number of charged tracks. The global cuts 
are quite loose in that they cause little (typically < 5%) signal rejection. The track 
cuts are more stringent and so the existence of one or more good high-do tracks is 
the most important property of a signal event. The relative importance of track 
and global cuts and a rationale for the general selection technique is discussed in 
Section lO. 



6.1 Cuts common to all selections 

Track cuts 

All selections require at least one good high-do track (N goo d > in the notation of 
Section |6.2| ). This section describes the track cuts that define what is meant by 
'good'. 

The momentum, d and Xbs (^ ne X 2 °f the track fit to the beam spot per 
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ee 


/i/i 


TT 


e/i 


er 


\IT 


p (GeV) 


low dj 
high dj 


19.0 
19.8 


20.0 
18.6 


4.7 
6.1 


14.9 
20.0 


10.0 
8.5 


8.1 
4.3 


d (cm) 


low dj 
high dj 


0.17 
0.50 


0.50 
0.50 


0.13 
0.50 


0.09 
0.50 


0.24 
0.49 


0.5 
0.48 


Xbs 


low 

high dj 


700 
680 


550 
690 


230 
700 


180 
700 


210 
590 


700 
610 



Table 6.1: The cuts on track momentum (p), do and Xbs f° r eac h channel under 
each optimisation. A track must have parameters greater than these values to be 
considered as a good high-do track in the corresponding selection. 



degree of freedom) are the three variables whose associated cuts are tuned, and so 
are required to be greater than values which depend both on the channel and on 
the optimisation in question. Table |6.1| shows these values. A track is allowed to 
fail the d cut if its z is > 8 cm. No other (global or track) cut depends on the 
optimisation, and so Table |6.1| completely summarises the difference between the 
low and high decay-length optimisations. 

The %i of a track is required to be greater than 1, and it must have at least 4 
assigned TPC hits and less than 3 bad TPC points (defined as a point where the 
track is missing a TPC hit for no observable reason - see Sections |5.5.1| and p. 5. 2 ). 

For the purpose of eliminating photon conversions the track has a variable called 
r po int- This is the radius at which the track's circular rip projection is pointing 
directly away from the origin. If the track has more than one hit at a radius below 
this then r po i nt is set to zero. Tracks originating from photon conversions should 
be pointing away from the origin at the point where the conversion occurred, and 
should have no hits at lower radii. Thus tracks with 185 cm > r point > 4.5 cm are 
rejected. This is almost never satisfied by signal tracks, and is a very effective cut 
against photon conversion tracks since it does not rely on both the electron and 
the positron from the conversion being reconstructed. 

To further reduce the background from nuclear interactions the distance of 
closest approach in three dimensions of each track to its nearest tagged nuclear 
interaction vertex (if there is one) is recorded and forms the variable d^j. Tracks 
with djsij < 5 cm are rejected. If the track is part of a multi-track vertex that 
failed both the nuclear interaction and f — * T mu i ti _ pron g criteria, then it is rejected 
if the x 2 °f the vertex is less than 100, and its maximum two-track angle is greater 
than 20° or its kink angle is less than 20°. The latter two conditions reduce the 
probability of rejecting vertices that are due to slepton decay but which have not 
been identified as such. 

To reduce the background from cosmic rays comparisons between the candidate 
tracks of an event are made. If two tracks do not have a reconstructed slepton- 
lepton vertex with a kink angle in rip greater than 90° at a radius greater than the 
4th TPC pad row, and the difference between their d 's is less than 10% of their 
average do, then they are both rejected. 
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Tracks with Olf = 5 are required to have a slepton-lepton vertex radius less than 
that of the 20th TPC hit pad and a slepton momentum greater than 2 GeV (to 
further reduce background from ECAL splash-backs). Also, the track-fit \ 2 (P er 
degree of freedom) must be less than 4 (this is the fit of the helical track to the 
hits, not of the two tracks to the TPC- vertex). 

Tighter cuts are made on tracks with < 5. They are required to have at least 
one VDET hit with r<p information or do> 2 cm or Zo> 10 cm. The track-fit x 2 is 
required to be less than 2.5. To reduce the cosmic ray background the momentum 
is required to be less then 150 GeV. To reduce the background from low angle 77 
events the | cos(0)| of the track is required to be less than 0.94, and if "Jij = 4 the 
I cos(6 l )| of the slepton-lepton vertex is also required to be less than 0.94. Cuts are 
also made on the minimum possible I — > I decay angle in two dimensions (the rip 
projection), a miri2 , and in three dimensions, a min3 . These are defined as the angle 
between the track and the line joining the track to the IP at the track's innermost 
hit location. It is required that a min2 > 1.5° and a min3 > 5°. 



Global cuts 

A certain degree of cosmic elimination is obtained by requiring that the event is 
not 'class 3'. Classes are groups of basic cuts that an event is tested against after 
reconstruction. There are 25 in total and all events satisfy at least one of them. 
They are intended to speed up analyses by providing fast rejection of events that 
a user wishes to ignore. Class 3 is intended to flag cosmic rays passing through the 
VDET. It requires exactly two tracks with d < 10 cm, z < 50 cm, | cos(6 l )| < 0.95, 
p > 10 GeV and with > 9 TPC hits. It then requires that at least one of these must 
have a do > 1 cm or a z > 2 cm, and that their directions differ by > 179° in 0, 
< 1° in 9, and their do's differ by < 0.5 cm. 

Further rejection of cosmic events is obtained by requiring to < 100 ns, and that 
either the Oij of the best track is greater than 4 or that the total charged momentum 
from the IP is greater than 1 GeV or the cosine of the angle between the two highest 
momentum tracks that have TPC hits is > —0.999. If the "Jij of the best track is 
less than 4 the additional constraint that the cosine of the angle between the two 
highest momentum tracks (regardless of hits) is also > —0.999. 

Beam-gas events are rejected by requiring that the charged momentum from 
the IP be non-zero or that the ratio of transverse to longitudinal momentum from 
the IP be > 0.4 (beam gas events are low angle and not centred in the detector). 



6.2 Additional cuts for individual selections 



In addition to the track cuts described in Section |6TT| , tracks are also subject to an 
additional cut if they are to pass the ee, jlfl or ep, selections. For the ee selection 
only tracks with a non-zero probability of being an electron are accepted. The 
probability is calculated during reconstruction using calorimeter information. It is 
equal to zero if the track passes some very loose muon identification criteria. For the 
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jljl selection only tracks with a muon probability (calculated during reconstruction 
using calorimeter and muon chamber information) greater than 0.2 are accepted. 
For the ep, selection a track must pass either of these conditions. No such cuts can 
be used in the channels involving staus since the high-do tracks can be electrons, 
muons or pions. 

Events are also subject to additional selection-dependent global cuts. The global 
variables that are used are described in Table |6.2| . Primed variables are calculated 
only using tracks and calorimeter deposits selected by the energy flow algorithm. 
For charged tracks this puts a cut of 2 cm on do and 10 cm on z since it is only 
intended for physics at the IP. The phrase 'from the IP' in Table |6.2| can be trans- 
lated as 'defined by energy flow'. Variables with an aTrow' above them pertain to 
the 'best' track in the event. This is defined as the track with the highest value of 
Olj. If there is more than one track with a of 5 (the maximum) then the best 
is the one with the highest momentum. All energies and momenta are in units of 
GeV, all angles are in degrees. 

Several variables refer to 'hemispheres' or 'jets'. The event hemispheres are 
formed by finding the thrust axis of the event. This is the unit vector, n, that 
minimises Yl P n ) where the sum is over all energy flow objects (charged tracks, 
photons and neutral hadrons) and p refers to the momentum 3-vector of a given 
object. Objects are then sorted into the hemispheres depending on which side of 
the axis they are closest to. The direction of a hemisphere is the direction of the 
sum of the 3-vectors of all its constituent objects. The jets are formed using the 
commonly used DURHAM algorithm [ fHfl . The y-cut is chosen such that the event 
has exactly two jets, unless the event has only one object in which case jets are not 
defined. The direction of a jet is the direction of the sum of the 3-vectors of all its 
constituent objects. 



Variable Description 

N c h Total number of reconstructed charged tracks 

N' ch Number of reconstructed charged tracks from the IP 

N e Number of electrons 

P e Probability that the best track is an electron 
Jlj The 3lj of the best track 

C2ch The cosine of the angle between the two highest momentum tracks 
Ngood The number of tracks passing the track cuts (must be > 1) 

S' 2j3 a/1 — 0.5((3f + /?|); where (3\ and (3 2 are the boosts of the two hemi- 
spheres (one of which can be zero if the hemisphere contains no energy 
from the IP) 

$' aco The acoplanarity of the jets (the angle between their rip projections; 

defaults to 999 if the event cannot be forced into 2 jets) 
Ntpc The number of tracks that have TPC hits (tracks that have been tagged 
as originating from a non-SUSY secondary vertex are not included) 

continued on next page 
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Variable Description 



N' 



a 



i 

hemi—z 



E tot 
~f 
P't 
P' Z 
2/23 



E, 



beam 



thrust 



E' 



12 



N, 



60 GeV 



m 
a 



i 

tot 
I 



tot 



Ptot 



The number of tracks from the IP that have TPC hits 
The minimum angle between the plane of the hemisphere directions and 
the z-axis; defaults to 999 if one or more hemisphere contains no energy 
from the IP 

The total energy from the IP 

The momentum of the best track 

The transverse momentum from the IP 

The longitudinal momentum from the IP (unsigned) 

The DURHAM algorithm y cut for which the event goes from 2 to 3 jets 

The | cos(0)| of the momentum missing from the IP 

The LEP beam energy (constantly updated for real data by the LEP 

Energy Working group using 15 minute chunks of data, see JOJ) 

The | cos(0)| of the thrust axis 

The I cos(60| of the best track 



The energy from the IP within 12 degrees of the beam axis (for Monte 
Carlo data this has an extra amount added taken randomly from a dis- 
tribution determined from random triggers in order to simulate beam 
noise) 

The number of tracks with do< 0.5 cm and zq< 3 cm and a dE/dx no 
less than 10 standard deviations below that which would be expected 
from a singly charged particle with a mass of 60 GeV (the approximate 
GMSB slepton mass limit) 
The invariant mass at the IP 

The acollinearity of the hemispheres (the angle between them; defaults 
to 999 if one or more hemisphere contains no energy from the IP) 
The number of muons 

The total charged energy (tracks that have been tagged as originating 
from a non-SUSY secondary vertex are not included) 
The total charged momentum (tracks that have been tagged as originat- 
ing from a non-SUSY secondary vertex are not included) 



Table 6.2: The variables used in the selections and their definitions. 



It may seem best to develop global cuts which, like the track cuts, tightly define 
the signal properties and so exclude background by the fact that it is non-signal-like. 
For example, it may seem sensible to require in the ee selection that 4 < N c h < 6 
(since there should be four lepton tracks with possibly two slepton tracks) and 
N e = 4. But such simple cuts, while intuitive, cause significant losses in the signal 
efficiency The reason for this is twofold. Firstly it is a matter of how large a loss of 
efficiency is considered significant. Because of the observable slepton lifetime the 
signal is extremely distinctive (see Figure |4.1|) , and so it should be possible to gain 
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close to 100% efficiency with little background as long as the decay-length is of the 
right size. Thus while a cut which reduces the efficiency by (for example) 5% may 
be considered perfectly acceptable in most SUSY searches, here it is considered to 
have a significant detrimental impact and alternative cuts would be sought. The 
other reason is that the signal is simply not well defined by global variables. The 
slepton lifetime makes the usually indispensable tool of energy flow of very limited 
value, and interferes with the particle identification since the lepton tracks can be 
foreshortened and can cross the TPC wires and enter the calorimeters at unusually 
high angles to the radial direction. To go back to the earlier example, the cut 
4 < N c h < 6 causes an average loss of efficiency in excess of 10% over all the ee 
Monte Carlo samples generated at 208 GeV, and just the cut N e > 2 causes more 
than a 10% loss in nearly half the samples. Thus to obtain the high efficiencies 
that must be achievable given the nature of the signal, the cuts on global variables 
must be set loose. 

The strategy then was to focus the global cuts more on the properties of the 
background rather than that of the signal, and to ensure that they had virtually 
no effect on the signal efficiency for the whole (m x , mj, dj) signal parameter space. 
Of course, in the case of slepton decay-length comparable to the do-resolution of 
the detector, some cuts will have to have a severe impact on the efficiency. But 
these cuts should be those related to the one, well-defined and principal feature 
of the signal - the slepton decay-length, and so should be track cuts. Also, if the 
number of cuts with the potential to harm the efficiency are kept low, they can 
all be tuned by an optimisation process and the analysis as a whole will then be 
close to optimal. These are the cuts on track momentum, d and Xbs^ which have 
already been described. 

So the most important cuts in the analysis (those cuts which reject the large 
majority of the background) are the track cuts, and the global cuts serve only as 
a means of further discrimination in the case that an event has at least one good 
high-do track. 

The additional global cuts are given in tabular form in Sections p.2.1\ to p.2.t . 
The third column in each table indicates which Standard Model background the cut 
is intended to remove. If a cut is quite general in that it is effective against many 
backgrounds, then a '-' is used. The cuts are grouped into numbered sets, where 
each set is intended to discriminate against the stated background. Ideally, every 
cut would have its associated plots of the distributions of the associated variable 
for signal and for background presented here, but given the number of cuts this is 
impractical, and so plots will only be shown for a few selected cuts. 



6.2.1 Additional cuts for the ee selection 

The cuts are given in Table |Q| . The event is required to pass cut 1, and is only 
required to pass cut 2 if N goo d = 1 and Jij = 2 (i.e. minimally signal-like - only one 
track passing the cuts with no evidence of the parent slepton). Cut 2 is an anti- 
rr cut that discriminates against the high energy and (frequently) back-to-back 
nature of rr events. It, or a slight variant of it, is used in every selection, and the 
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distributions of the associated variables for the ee channel and rr background are 



shown in Figure |6.1| on page |83]. rr is a major background due to the r lifetime 
which gives it a decay length of ~ 5 mm at Ef, eam ~ 100 GeV, and the fact that the 
multi-prong decays are narrow which can lead to reconstruction difficulties as hits 
can be incorrectly assigned amongst the tracks from a single decay. 
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20>N ch >2 


K h <7 


{iV e >0 or P e > 0.1999 


or %>4} 


qq 
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0.99>c 2ch > 


-0.999 


{^>0.1 or 


$' < 174 j 

aco J 


TT 



Table 6.3: Additional cuts for the ee selection 



6.2.2 Additional cuts for the jljl selection 

The cuts are given in Table |6.4| . The event is required to pass cut 1. As with the 
ee selection it is only required to pass cut 2 if N good = 1 and = 2. 
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20>N ch >2 




Kh<7 
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0.99 >c 2ch > -0.999 


{^>0.1 


or $' aco <174} 
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Table 6.4: Additional cuts for the fifi selection 



6.2.3 Additional cuts for the ff selection 

The ff selection is, as might be expected, the most complex of the six. The high-do 
tracks can be electrons, muons or pions and so no additional cuts can be made on 
the track properties. There are also difficulties with the global properties of the 
events. Whereas in a channel without staus there are only two invisible particles, in 
a ff event the neutrinos from tau decay mean that there can be up to ten, leading 
to a much greater spread in the visible energy, and increased background from 77 
events. The multi-prong decays also push up the multiplicity of the events which 
make qq events more difficult to remove. 

The cuts are given in Table |6.5| . The event is required to fail cut 1. This is 
because cut 1 is intended to pick out 77 events from f f events rather than the other 
way around. If % = 5 then this is the only cut the event has to pass, although 
this is very rare in background events. If 3lj < 5 then the event also has to pass 
cuts 2 and 3 which are simple cuts on energy, mass, multiplicity and acollinearity. 
If there is no f — > T mu iti- prong candidate in the event it is also required to pass 
cuts 4 and 5, and to fail cut 6 (another pro-77 cut). Under the hypothesis that 
the event is 77 and that only the electron/positron travelling in the direction of 
the missing momentum was deflected, then p' t / (E beam + p' z ) is the tan(#) of that 
electron/positron after deflection. Finally if % < 5 and there is no f — > T mu i ti _ prong 
candidate and o>' hemi _ z < 6 and N good = 1 and Olj < 3 then the event is also 
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required to pass cut 7, which is a variant of the anti-rr cut used in the previous 
two selections. 

Cut 4 is the most convoluted of the cuts. This is because it uses variables whose 
distributions for signal and background overlap significantly, but since Ntpc and 
a 'hemi-z are not correlated with each other or with c' thrust and c' miss (although c' thrust 
and c' miss are correlated), requiring events to pass only one of the cuts on these 
variables means the overall cut has significant impact on qq events (approximately 
90% fail) while leaving signal quite untouched (approximately 1% loss in efficiency). 



Figure [J72] on page [84] shows the distributions of these variables and the cut values. 
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qq 
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Table 6.5: Additional cuts for the rr selection 



6.2.4 Additional cuts for the efb selection 



The cuts are given in Table |6J| The event must pass cut 1 and, if N goo d = 1 and 
ftj = 2, cut 2 also. 
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20>iV c/i >2 N' ch <7 
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Table 6.6: Additional cuts for the e\x selection 



6.2.5 Additional cuts for the ef selection 



The cuts are given in Table |677]. The event is required to pass cuts 1 and 2. Cut 
3 is only applied if 5ij < 3 since it can be damaging to signal efficiency when the 
slepton decay length is long. In particular, an upper limit on E' tot is not applied 
if it looks like there is a fully reconstructed slepton in the detector (iV 60 GeV > 
and Ntpc < H) since the slepton's energy will be included by energy flow and so 
there may be no missing energy. Cut 4 is the standard rr cut and is only applied 
if Ngood = 1 and ftj = 2. 
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Table 6.7: Additional cuts for the er selection 



6.2.6 Additional cuts for the jif selection 



The cuts are given in Table p78| . The event is required to pass cuts 1 and 2. Cuts 
3 and 4 are the same as in the ef selection and are only applied under the same 
circumstances. Cut 5 is an additional anti-77 cut and is only applied if N goo d = 1. 
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Table 6.8: Additional cuts for the \xf selection 



6.3 A note on the exclusivity of the selections 

It should be noted that the selections are not mutually exclusive. That is to say, 
a single event can pass more than one of them. In fact, there is no combination of 
selections which an event, in principle, cannot satisfy. This is because the channels 
all give very similar final states, and so reliably disentangling one from another is 
practically impossible. 

But exclusivity is very desirable amongst the selections of a search with more 
than one channel, since the results of the selections are then statistically indepen- 
dent, making any subsequent combination far simpler. Deliberate efforts could be 
made to ensure the selections are exclusive, but the method would have to be quite 
arbitrary, and would have a similarly arbitrary effect on the discovery/exclusion 
power of the search. 

It would of course be possible simply to have a single selection aimed at selecting 
signal events from any channel. But in order to encompass all six channels without 
reproducing the same level of complexity of six individual selections, an all-inclusive 
selection would necessarily have to have a higher level of background acceptance. 

The final states of the channels of the signal process fall mid-way between 
being different enough to allow exclusive selections, and being similar enough to 
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allow a single selection, without harming the final result. Individual non-exclusive 
selections have been used then, and the complexity that follows is dealt with in the 
next chapter. 

6.4 Optimisation 

As has already been stated, the cuts on track momentum, do and Xbs are the mos t 
critical of all the selection cuts. They cannot be set loose to encompass all the signal 
for any decay-length since they are very important in rejecting background, but 
if set too tight they will have a severe detrimental impact on the signal efficiency. 
They are therefore tuned to obtain the optimal cut values. 

The definition of the optimal cut value is that which gives the lowest average 
expected 95% confidence level upper limit on the signal cross-section, a 95 , that is 
the lowest average 95% confidence level cross-section limit that would be obtained 
over a large number of experiments in which the outcomes were as expected under 
the background-only hypothesis, as dictated by background Monte Carlo studies 
(the real data is not used). The methods for calculating the confidence level and 
<7 9 5 will be given in Sections |0| and [T]5] respectively. 

Ideally a 95 should be obtained as a function of all the cut values simultaneously. 
But since there are eighteen of these (three variables x six channels), and the 
method for calculating just a single cross-section limit is quite complex and CPU- 
intensive, this is impractical. Thus (J95 is obtained as a function of each cut value 
individually, first in the selection order ee, jip,, ff , ep,, ef, p,f for the cut on 
momentum, then the same for the cut on Xbsi anc ^ t nen f° r the cut on do- At each 
stage the cut value that gave the smallest 095 was identified and the cut set there. 
Cuts waiting to be tuned were set at sensible initial values. 

A (X95 can only be calculated though, for a given signal hypothesis. That is, 
a given point in the (m x , mj, df) parameter space and a given set of values for 
the x — > e, /i, f branching ratios (Bg^?). Eight such signal hypotheses were used, 
corresponding to all permutations of high and low decay-length, high and low 
neutralino mass, and either slepton co-NLSP (-Bg,/i,f = §) or stau-NLSP {B^ = 0, 
Bf = 1) branching ratios. The high and low decay-lengths used were 80 cm and 
0.2 cm, and the high and low neutralino masses were 95 GeV and 88 GeV. 

It was found however, that good sensitivity could be obtained under all the 
signal hypotheses using just the two optimisations obtained from the low decay- 
length, low neutralino mass, slepton co-NLSP combination, and the high-decay 
length, high neutralino mass, stau-NLSP combination. Since it was found that 
the optimised cut values were most strongly influenced by the choice of slepton 
decay-length, these optimisations were dubbed the low and high decay-length op- 
timisations respectively, and the other optimisations were not used. 
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Figure 6.1: The distributions of c 2c h, S' 2f3 and §' aco for rr and ee Monte Carlo events 
at 189 GeV. Top left is the full distribution of C2 C h fc> r TT events. The two other 
plots on the same row show sub-ranges of the full distribution and the cut positions 
used in the anti-rr cuts. The second and third rows show the same for ee events 
with decay-lengths of 0.2 and 500 cm respectively. The fourth row shows scatter 
plots of $' aco versus S' 2 p for the same samples. The rectangle in the top left of each 
plot shows the anti-rr cut on these variables (events inside the rectangle fail the 
cut). 
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6.2: The distributions of E' 12 , Ntpc, 



a. 



and of c'.i 



versus c: 



for 



hemi—z aLLK ~ L KJL '"thrust ••"'-^•^■^ ■"miss 

ff Monte Carlo events at 189 GeV. These are the variables used in cut 4 of 
selection. In the distributions of a' hemi _ z , values in excess of 10 have been 
in the last bin to retain the overall normalisation. The plots on the left 
qq events, the centre plots for ff events with a decay-length of 0.2 cm, and 
ht plots for ff events with a decay-length of 500 cm. The cut values used 
ff selection are shown as lines in the top three rows of histograms, and as 
des in the bottom row (events outside the rectangle fail the cut). 



Chapter 7 

Results and interpretation 



7.1 Outcome 



The number of background events expected to pass at least one selection, calculated 
from applying the selections to the background Monte Carlo samples, is shown in 
Table |7.1| for both optimisations (see previous chapter) along with the number of 
events observed in real data. The data is consistent with the background expecta- 
tion, and so there is no evidence of supersymmetry. The numbers are broken down 



into seven bins by LEP energy - a re-binning with respect to Table Since Monte 
Carlo background studies were only performed at energies of 189 and 208 GeV, the 
expected numbers of background events for intermediate energies were gained by 
re-normalising the values at these two energies to the correct luminosity, and then 
linearly interpolating between them. 
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206.7 


133.7 


1.2 ±0.5 


2 


0.3 ±0.2 


1 


Total 


5.2 ±i.4 


5 


1.5 ±0.8 


4 



Table 7.1: A comparison of the expected number of events from background Monte 
Carlo studies with the number observed in data, broken down by LEP energy. L is 
the total integrated luminosity taken at each energy. The quoted errors are purely 
statistical in origin. 



Since the selections for each channel are not mutually exclusive, a channel-by- 
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channel breakdown of the number expected versus observed is not appropriate. 
Rather, the numbers (expected and observed) are broken down into those passing 
particular sets of selections. Since there are six selections there are (2 6 =) 64 
possible outcomes for the combination of selections that any event passes, including 
the complete set (passing all selections) and the null set (failing all selections, and 
therefore of no interest). The selection sets are numbered from to 63 according 
to the convention 

So-^63 = 32 5gg + 16 5^ + 8 Sff + 4 5ejl + 2 Sgf + Sfif , 

where 8 c h<m is for events failing the respective channel's selection, and 1 for 
events that pass. The expected versus observed number of events is shown broken 
down by selection set in Table |7]2|. Only sets for which the expected background 
was non-zero, or for which a data event was observed are listed. No errors are 
quoted on the expected numbers of events in this table since, by breaking the 
numbers down by selection set, the events passing the selections become so thinly 
spread that each figure is typically based on no more than one event; and the 
procedure for calculating the (asymmetric) error on a weighted sum of Poisson 
variable estimators, when the estimators themselves are based on very few events, 
is extremely complex. Even without errors though, it is clear that the agreement 
for sets with a non-zero expected number of background events is good. Two events 
under each optimisation are observed however, with sets under which there were 
no events expected. Given that, overall, the number of background Monte Carlo 
events analysed was only around three times the number of events expected in 
the data set however, this does not constitute an incompatibility between real and 
background Monte Carlo data. 



7.2 Efficiencies 



Figures [Fj] to [7.31 show, for each channel, the probability of a signal event passing 
its own channel's selection, and the probability of it passing at least one selection 
no matter which, as a function of decay length (estimated using the prescription 



described in Section [4.2. 1| ) for both the low and high decay length optimisations. 



All plots are created from the Monte Carlo data generated at y/s = 189 GeV. 

It can be seen that, in general, the selections other than that of the channel's 
own selection add little to the total efficiency for a channel (i.e. if an event is not 
selected by its own channel's selection, it is unlikely to be selected by another), as 
might be expected. 

In each plot the different efficiency values at each decay length correspond to the 
different points in (m x , mj) space at which Monte Carlo data was generated (listed 
in Appendix [A]). All the points associated to the same {m x , mj) point are joined 
by a dotted line. This is for guidance only, it does not represent an interpolation. 
The dashed lines are formed by the straight line extrapolations from the two lowest 



and two highest decay length points. More will be said about these in Section O. 
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0.0 
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14 


0.0 
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Table 7.2: A comparison of the number of expected events from background Monte 
Carlo studies with the number observed in data, broken down by selection set. 



Note that there are many more mass points per decay length for the ee channel 
since it was used as a test case to determine how many points in mass space would 
be needed for the subsequent Monte Carlo samples. 

The very weak explicit dependence of the efficiencies on the neutralino and 
slepton mass is manifest. Most of the variation in the efficiency between the mass 
points at a given decay length can be accounted for by statistical fluctuation from 
a common value. In each plot there is one obvious exception however, most clearly 
at high decay length. In every plot the point with the highest efficiency at the 
highest decay length is the point with the least initial (i.e. before slepton decay) 
phase space open to the events. Specifically, for channels not involving staus it is 
the (m x , mj) point (94GeV, 93GeV); for channels that do involve a stau it is the 
point (94GeV, 91GeV). In both cases this is only a 0.5 GeV difference between 
the beam energy and the neutralino mass, and approximately a 1 GeV difference 
between the neutralino mass and the mass of the slepton plus that of the (heavier) 
lepton. The reason for the larger efficiency at high decay length for these points is 
that their decay length is greatly overestimated. The lack of phase space results in 
the slepton moving very slowly, so that its dE/dx is very large. Thus the slepton 
is stopped by energy loss through ionisation almost as soon as it encounters any 
matter, and so typically does not penetrate the beam-pipe. Thus its decay length 
is effectively limited to ~ 5 cm. The efficiency does peak and then drop with decay 
length in a not too dissimilar way to the other points, but the drop-off at high decay 
length is not because the sleptons are increasingly escaping the tracking volume 
before decay. It is because in order to give such a slow-moving slepton such a 
large (estimated) decay length, the slepton lifetime has to be set very large (0.2 /is 
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for a selectron at the stated mass point with dj = 500 cm). Once it is stopped 
in the beam-pipe it waits out the remainder of its lifetime before decaying. This 
can exceed the time period during which the ITC is active (~ 1/xs), and so the 
resulting decay products (the high-do tracks) do not have ITC hits and so fail the 
quality criteria described in the previous chapters. This is increasingly likely with 
increasing estimated decay length and so causes the associated drop in efficiency, 
but still results in better efficiencies than the other points for the highest decay 
length. 

The efficiencies for the = 208 GeV Monte Carlo data are not given since 
they are virtually identical when compared with the appropriate scaling (more will 
be said about this in Section |7.4j ). The slow-slepton effect resulting in the higher 
efficiencies for high decay lengths is not observed in the 208 GeV Monte Carlo data 
though, since the minimum E beam — m x energy difference was limited to 1 GeV. 
This shows that the effect is only confined to a very small corner of the mass space. 

7.3 The confidence level 

The ultimate goal of a search, in the absence of observing significant evidence 
for the signal, is to place the most stringent possible limits on the parameters of 
the underlying theory. Due to the statistical nature of particle physics, a search 
can never, with 100% certainty, rule out the hypothesised phenomenon. Only a 
'confidence level', necessarily less than 100%, can be attached to the assertion that 
the phenomenon was not occurring. The confidence level is the probability that, 
if the hypothesised process were actually occurring, and if the experiment were 
done again in exactly the same way, a result more compatible with the hypothesis 
would be obtained. A high confidence level then indicates an experiment that is 
not compatible with the hypothesis, a low confidence level indicates one which is. 
The benchmark confidence level for 'exclusion' (disproof of the hypothesis) is 95%. 
At each point in the relevant parameter space the confidence level can be computed 
from the slepton and neutralino properties, and those points for which it is > 95% 
are said to be excluded. This section describes the way in which the confidence 
level is constructed. However, the method of construction used for this analysis is 
unusually complex for a SUSY search, and so there is first a discussion highlighting 
the motivation for this complexity. 

The confidence level is computed from a 'test statistic', t, which itself is com- 
puted from the information yielded by the search. The test statistic should be 
a measure of how 'signal-like' (how compatible with the search hypothesis) the 
search outcome is. The higher its value, the more signal-like. Then if the expected 
probability distribution function of the test statistic under the signal hypothesis, 
p(t), is known, the confidence level is given by 
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Figure 7.1: Efficiencies for the ee and \x\x channels as a function of estimated slepton 
decay length. The different curves correspond to different points in the (m x , m\) 
mass space. 
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Figure 7.2: Efficiencies for the ff and eji channels as a function of estimated slepton 
decay length. The different curves correspond to different points in the (m x , mj) 
mass space. 



7.3. THE CONFIDENCE LEVEL 



91 



_i5 0.9 
'o 

HE °- 8 

LJ 

0.7 

0. 

0. 

0. 

0. 

0. 



ex 



ex 



ex selection I 
for low decay-length 



1 -o 10 10 

Estimated slepton decay length (cm) 



Any selection !:% 
for low decay-length \\\ 

1 10 1 o 2 1 o 3 

Estimated slepton decay length (cm) 



0.6 
0.5 
OA 
0.3 
0.2 



: 0.2 0.5 2 


■ 


300 


ex 










j 


ex selection 








I for high decay-length 

















ex 



Any selection 
for high decay-length 



10 1 '0 10 

Estimated slepton decay lengt 



1 10 10 10 

Estimated slepton decay length (cm) 



0.2 0.5 2 5 



0.2 0.5 2 5 



50 80 300 



|IX 



[IT 



\xt selection 
for low decay-length 



1 10 1 o 2 1 o 3 10 

Estimated slepton decay length (cm) 



Ally selection 
for low decay-length 



1 10 10 10 

Estimated slepton decay length (cm) 



fix 



ixx selection 
for high decay-length 



o.i 





Any selection 
for high decay-length 



10 1 10 10 10 

Estimated slepton decay length (cm) 



1 10 10 10 

Estimated slepton decay length (cm) 
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mass space. 
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where t 



■mm 



is the minimum possible value of the test statistic, and t 



meas 



is the 



actual value from the search result. The definition of t is, formally, completely 
free. Clearly it should be constructed such as to have good discriminating power 
between the signal and no-signal scenarios, but a bad choice only leads to a non- 
optimal confidence level (lower than the best possible). It does not make the 
confidence level wrong. 

For the most optimal confidence level, p(t) should take into account that some 
contribution to the value of t is expected from background as well as from signal, 
and so should be the density of the test statistic under the signal-plus-background 
hypothesis. Any systematic errors in the simulation of the background however, 
then find their way into the confidence level. So the background contribution to p(t) 
should be left out unless the background is well understood and its simulation is 
trusted. Since the background to this search comes from unconventional processes 
(principally particle interaction with the structure of the detector), the simulation 
of which have not been checked with the same rigour as more conventional sources, 
the background contribution to p(t) is left out. 

Typically with searches, as is the case with this search, the larger the number 
of observed events, the more signal-like the outcome. If there is no additional 
information associated with each event to discriminate signal from background, 
the only sensible choice for the test statistic is then just the number of observed 
events, n obs . The confidence level is then given by 



where P n is the (Poisson) probability of obtaining n events as a search outcome, and 
n exp is the mean number of expected events under the signal (with no background) 
hypothesis. 

In this way, a confidence level can be obtained for each channel by using 



where L is the integrated luminosity, a is the neutralino cross-section, P c han is the 
probability that the neutralinos will decay such as to give this channel (given by 
the branching ratios), and e c h an is the efficiency with which events of this channel 
will be selected. For statistical reasons it cannot be said that the point is excluded 
if any one of these CLs is greater than 95%. The channel that appears most likely 
to exclude the point based on Monte Carlo studies only should be chosen, and if the 
CL of that channel (and that channel alone) is greater than 95% then the relevant 
point in parameter space is excluded. 

This is clearly, however, a non-optimal approach. The information from five of 
the six channels is being thrown away. It should be possible, rather than to perform 
CL calculations for each channel individually, to calculate just one CL using all the 
data which would be greater than any one of the single-channel CLs. This cannot 




n=0 



n=0 



n bs = number of events passing channel's selection, 

f^exp & @ Pchan ^chan j 
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be calculated in a simple way from the single-channel CLs directly because the fact 
that the associated selections are not mutually exclusive means that the CLs are 
not statistically independent. A solution would be to use 

n b s = number of events passing any selection, 

chan=l 

But this is still non-optimal because it treats all observed events in an equal 
way. The number of expected signal events will be different between channels, and 
the number of expected background events will be different between the selections. 
Thus the contribution of an event to the value of the test statistic should be in 
some way representative of the discriminating power of the selection(s) it passes. 
In particular, in the case that the stau is significantly lighter than the selectron and 
smuon, and so effectively becomes the sole NLSP, it can be expected that including 
events that do not pass the ff selection will only worsen the result. 

A solution is to use the non-exclusivity of the selections as an extra source of 
information. The distribution of the signal efficiency amongst the 63 selection sets 
(disregarding the null set) is quite different from the equivalent distribution of the 
expected background. Thus the set of selections an event passes can be considered 
a discriminating variable that gives information on how signal-like it is. If the sets 
are numbered we have an integer variable, denoted S, defined for each event with 
a value between 1 and 63. The test statistic can then be defined as the ratio of 
the probability of the search outcome under the signal-plus-background hypothesis 
to that of just the background hypothesis (i.e. a relative likelihood) taking into 
account the known numbers of signal and background events expected with each 
selection set: 

n^( s s + &s,^s) 
_L f± 6 =J 

s 

P(x, n) is the Poisson probability of observing n events from an expected mean 
of x, s s and b s are the mean numbers of signal and background events expected 
to pass the set of selections, S, respectively, and n s is the observed number of 
events passing selection set S. This test statistic now makes use of all the available 
information and copes naturally with variations in the number of expected signal 
events per channel and number of expected background events per selection. 

A drawback is that p{t) cannot be determined analytically now, only by a 'toy' 
Monte Carlo experiment. That is to say, using random numbers, a succession of 
hypothetical search results has to be generated in accordance with that which would 
be expected under the signal-only hypothesis. The corresponding test statistic can 
be calculated for each and a normalised histogram of these values built up to form 

tit)- 

The sorting of events by selection set now requires an associated decomposition 
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of the efficiency. Thus we refer to the 'total' efficiency, defined as the probability 
of a signal event passing at least one selection, and the 'partial' efficiencies, defined 
as the 63 probabilities of a signal event passing each respective selection set. Then 

Sg = £j(J ^ Pchan^S chart > 
chan 

and 

63 

£-chan ^ ^ chan 
S=l 

where e s denotes the partial efficiency associated with selection set S, and e denotes 
the total efficiency. 

So far the discussion has not incorporated the fact that the search detailed in 
this thesis spans a range of energies at which different amounts of data were taken, 
and for which the signal cross-sections and partial efficiencies will be different. This 
should also be built into the confidence level. If the data is split into a number of 
energy bins then the total test statistic is just the product of the single-energy test 
statistics calculated for each bin. Thus the generalisation of Eq. [7.1| to multiple 
energies is 

'=n^nn \y . <-> 

where e denotes the energy bin and 

S Se = L e a e J2Pchane S chanc 
chan 

All the terms have the same basic definitions as before, but now a subscript e means 
that value as defined at the centre-of-mass energy corresponding to energy bin e. 

It is important to note that although numbers based on the background ex- 
pectation are used to calculate the confidence level, any systematic error in these 
numbers can only, on average, reduce the confidence level. This goes back to a 
point mentioned at the start of this section - that the definition of the test statistic 
is completely free. A bad choice for the test statistic (or using incorrect numbers in 
the calculation of it) merely yields a less than optimal number. It is only in the cal- 
culation of the density of the test statistic under the signal hypothesis, p(t), that 
systematic errors can cause a systematically inflated value, and the background 
expectation is not included in this. 

For a more detailed discussion on the construction of test statistics, including 



the case of possible systematic errors in the background expectation, see [[4(| . 
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7.4 Interpolation and extrapolation 

This section describes the methods of interpolation and extrapolation employed to 
obtain the signal efficiencies for any given point in (m x , mj, dj, y/s) space. Such 
methods are required for the scans that will be described in Sections |7]5| and |7.6| . 
The process is performed for each channel individually in stages. It is the partial 
efficiencies that are interpolated, and the total interpolated efficiency is their sum. 
Interpolation is first performed for four points in mass space to the correct mass 
point at decay lengths and y/s's above and below the desired values, then in decay 
length to form two points with the correct dj, and then in y/s to arrive at the 
desired point. The following describes the procedure for a single channel. 



Interpolation in mass space 

For a given slepton decay length and y/s, a method to gain the efficiencies for any 
point in (m x ,m ; -) space was required. It was decided that a linear method was 
preferable for reasons of simplicity and transparency. The simplest method for 
linear interpolation in two dimensions is bilinear interpolation, but since the points 
in mass space do not form a rectangular grid, bilinear interpolation would lead to 
discontinuities in the (partial and total) efficiencies as a function of the masses. To 
get a continuous function the mass points have to be joined to make a tessellating 
set of polygons, and then an interpolation performed over each polygon. Since 
there is no obvious benefit from using polygons with a larger than minimal number 
of sides, only added complexity and less transparency, triangles were chosen as the 
polygons into which to split the plane. A commonly used and natural method 
for dividing a plane into triangles based on a set of points is called Delaunay 
triangulation (originally described in |S7|), and this was the method used. 



The Delaunay triangulation of a point-set is the collection of triangles which 
satisfy an "empty circle" property, i.e. the circle drawn through the vertices of any 
triangle does not enclose any other point from the set. As long as no circle can 
be drawn through four or more points then the collection is unique. If this is not 
the case then there are degenerate collections (i.e. more than one collection that 
satisfies the Delaunay criteria). This is the case with the sets of points in (m x , mj) 
space considered here, which do have subsets of four points (but never more) lying 
on common circles. Thus the degeneracy necessitates a choice for the diagonal down 
which the quadrilaterals defined by the co-circular points are to be split in order to 
form the triangles. The decision is always taken to split such quadrilaterals down 
the diagonal with the lowest average total efficiency in order to be conservative. 
The triangulation is thus performed for each channel, energy, decay length, and 
optimisation separately. It should be noted that no extrapolation is done in mass 
space - the efficiency is assumed to be zero outside the convex hull of the point-set. 



Figure |7.4| shows examples of interpolated surfaces created by this method, 
and their associated Delaunay triangulations. Note that the range of the vertical 
scale has been limited to accentuate the variation of the surface. If viewed with a 
vertical scale from to 1 both the surfaces appear almost flat. It is the fact that 
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the efficiency is in general such a weak explicit function of the masses for a given 
decay length that allows the interpolation in mass space to be performed with so 
few points. Exceptions exist however, at high decay length. It can be seen from 
the efficiency plots of Section |7]2|, that when the slepton mass is sufficiently close 
to the neutralino mass, and the neutralino mass is sufficiently close to the beam 
energy, the efficiency fails to fall as quickly with high and increasing decay length 
as for the other points in mass space. Thus the efficiency takes a sudden rise in 
the relevant corner of the mass space at the largest decay length. Since few points 
exist to map this rise an interpolation in this region would be dubious. In order to 
deal with this, the point with largest total efficiency at the maximum decay length 
has its total efficiency artificially reduced to that of the point with the second 
highest total efficiency, and the partial efficiencies are scaled accordingly. This was 
judged preferable to the large amount of extra Monte Carlo generation that would 
be needed to map this rise in such a small corner of the parameter space. 




Figure 7.4: The total efficiency as a function of neutralino mass and slepton mass 
for the ee channel at 208 GeV under the high-c^- optimisation at decay lengths of 
a). 30 cm and b). 500 cm. The highlighted points on the surface show the points at 
which Monte Carlo data was generated. All other points on the surface are inter- 
polated from these using the Delaunay triangulation which is shown superimposed. 



Interpolation in decay length 

As can be seen from the efficiency versus decay length plots of Section |7T^, a good 
method of interpolation in decay length would simply be to join the associated 
points with straight lines, and this is what is done. Although this is linear interpo- 
lation, it is linear with respect to log(d z ~) rather than d\ itself. Unlike the treatment 
of mass space, extrapolation is used to extend the region of sensitivity to decay 
lengths beyond the maximum and minimum values at which Monte Carlo data was 
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generated. This is performed by simply extending the straight lines joining the two 
highest and two lowest points in decay length down to the X - clXlS, clS shown in the 
plots. 

Extrapolation is possible because it is a mathematical certainty that the effi- 
ciency cannot drop completely to zero no matter how high or low the decay length. 
The efficiency can only approach zero asymptotically as the decay length falls to 
zero or extends to infinity. Thus the roughly straight lines that form the left and 
right sides of the efficiency 'humps' must level out at more extreme values. The 
extrapolations will thus give estimated efficiencies lower than the true efficiencies, 
and the extension of the lifetime range will be conservative. The beginning of 
this curve on the low decay length side can be seen in the plots for the high-d^ 
optimisation. 

There is an obvious problem with one (m x , mj) point in the /if channel under 
the low decay length optimisation however, since the extrapolated line at high decay 
length does not descend to the x-axis. For one plot it is exactly horizontal, for the 
other it actually increases. But this is fixed by the reduction of the efficiency of 
this mass point at its highest decay length that is performed in order to make the 
interpolation in the mass plane valid (as described in the previous section). The 
plots of Section [T]2] are shown without this reduction. With the reduction the line 
drops to zero with the other mass points. 



Interpolation in y^i 

The method of interpolation in y/s (the LEP centre of mass energy) is again linear, 
and is based on an empirical observation. The mass points chosen for the generation 
at 208 GeV were, for the most part, simply scaled from those chosen at 189 GeV. 
That is to say, while the points occupy different positions in (m x , m z ~) space, they are 
at approximately the same positions in ( F m * , ^-) space (where E beam = J~s/2). 
It was observed that, for any given channel and decay length, the total and partial 
efficiencies of a point generated at 189 GeV were virtually identical to those of 
the point generated at 208 GeV which was closest in ( F m * , — ) space with the 
same decay length. This simple scaling indicated that no further Monte Carlo data 
was required at intermediate energies - the efficiencies could be obtained for any 
intermediate energy by an appropriate interpolation. 

To obtain the efficiencies at a given decay length d~ { and a given mass point 
(m x , mj) at LEP energy Eb ea m, the equivalent scaled masses (m x , m~) at 189 GeV 
and 208 GeV were found using the prescription: 

m x = m 'x 

Ebeam ^beam 

and 

m x m' x 
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where E' beam is either 189 GeV or 208 GeV to obtain the respective equivalent 
masses, and mi ept on is the average of the masses of the two leptons involved in 
the channel (e.g. (m M + m T )/2 for the jlf channel). The inclusion of the lepton 
masses in the scaling protects the m x — m ; ~ mass difference from dropping below 
the mass of either of the leptons when scaling down in energy. Then the efficiencies 
can be obtained at the scaled mass points for the relevant decay length for both 
189 and 208 GeV, and then linearly interpolated to the correct energy, although 
the interpolation is not very important since the values will be very similar. 

A method of quantifying the agreement of the partial and total efficiencies 
between 189 and 208 GeV was developed. For each point in (m x , mj, dj) space 
at which 208 GeV Monte Carlo data was generated, the 189 GeV point closest in 
the ratio-space defined above with the same decay length was found. Note that 
although 360 points in total were generated at 208 GeV, only 333 (> 90%) are 
directly comparable to 189 GeV points because the same decay lengths were not 
always generated at both energies. For each pair two x 2 ' s were calculated, one as 
a measure of the agreement of the total efficiencies and one as a measure of the 
agreement of the partial efficiencies without reference to their sum (it is thus a 
measure only of the relative sizes of the partial efficiencies, with any disagreement 
of the absolute sizes caused by a disagreement of the total efficiency removed). The 
former is defined as: 



£j is a total efficiency, N is the number of events out of the total sample of 500 
passing at least one selection, and the subscript i refers to either 189 or 208 GeV 
Monte Carlo data depending on its value. The x 2 from the partial efficiencies was 
defined in a similar way for each selection set S, and summed over, but now since the 
total number of events being considered is not fixed (events failing all selections, 
500 — Ni of them, are ignored) the errors are treated as Poisson as opposed to 
binomial: 




where 



£iV 2 + £2^1 



V 1 + V 2 



500 ' 





63 



2 



where 



*7s = 



?7siVs2 + VszVsi 
V S i + V s2 




esi is a partial efficiency corresponding to selection set S, N S i is the number of 
events out of the total sample passing selection set S, and the other terms have 
the same definitions as for the previous \ 2 ■ M, the number of degrees of freedom 
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to which the \ 2 is normalised, is taken as the number of selection sets for which 
Nsi + N s2 > 0, as opposed to the total of 63. This is because many of the selection 
sets have zero associated events at both energies and so do not contribute to the \ 2 - 
If they are included in the number of degrees of freedom the x 2 ends up artificially 
low. 

Both of these x 2 's are plotted in Figure for the low and high dr optimisations. 
Good agreement between the relative sizes of the partial efficiencies is observed for 
all pairs of points. Good agreement between the total efficiencies is observed for 
most pairs of points, but some have x 2 's well in excess of 10 (seven are in excess of 
100), and these have been placed in the last bins of the respective plots. All these 
pairs contain the 189 GeV low phase-space 'problem' points described in Section [772. 
It is not unsurprising that disagreement is observed here, since the same level of 
Ebeam — ™x degeneracy is not replicated in the 208 GeV Monte Carlo data, and 
the comparison is performed without the previously described efficiency reduction. 
The difference in the total efficiencies can thus be taken as a result of their slightly 
different positions in the mass-ratio space in a region where an anomalous effect is 
causing rapid change. 

In summary, the signal efficiencies are almost imperceptibly weak functions of 
y/s for a given point in the stated mass-ratio space and a given dj. Thus effi- 
ciencies can be reliably gained at all energies between 189 and 208 GeV by linear 
interpolation of the efficiencies calculated at these two energies. 



7.5 Model independent cross-section limits 

The cross-section limit is defined here as the 95% confidence level upper limit on the 
lightest neutralino production cross-section at a specific centre of mass energy for 
a given set of the relevant sparticle parameters (the neutralino mass and branching 
ratios, and slepton masses and lifetimes). Cross-section limits cannot be calculated 
analytically. A series of trial values are chosen and the confidence level is computed 
for each. Successive values are picked so as to converge on a computed confidence 
level of 95%. Once a confidence level sufficiently close to 95% is attained, the trial 
value can be taken as the cross-section limit. 

Since the cross-section (a) is a function of y/s however, and the search spans a 
range in y/s, a different cross-section value must be used at each. These are related 
by the function, a (y/s), which defines the evolution of the cross-section with centre 
of mass energy. But the use of a specific model to gain cr(yfs) will lead to limits that 
are not universally valid. This problem can be circumvented if the cross-section 
limit is calculated at (or above) the highest energy considered in the analysis, 
y/s max , and the cross-sections at lower y/s are obtained by a linear interpolation 
of the cross-section at the highest energy down to zero at threshold. Then the 
cross-section limits are valid as long as a = and d 2 a/d(y/s) 2 is negative at 
y/s = 2 m x and there is no more than one stationary point in <j(y/s) for y/s < y/s max 



see Figure [7.6|) . Then if the real cross-section is larger than the cross-section 



limit at y/s max , it is also larger at all lower energies. This underestimation of 
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Figure 7.5: The x 2 distributions resulting from compatibility tests between 189 and 
208 GeV Monte Carlo data. In the two lower plots not all pairs of points gave a x 2 
in the range shown. Such pairs have been placed in the last bin. Good agreement 
is shown amongst all the partial efficiencies, and most of the total efficiencies. 
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the cross-section at lower energies for a given cross-section at y^s = \fs max leads 
to an underestimation in the confidence level, and thus an overestimation of the 
magnitude of the cross-section at \fs max needed to give a confidence level of 95%. 
Thus the cross-section limit, quoted at ^/s = \fs max , will be conservative under 
these basic assumptions. 



Cross-section, o 

A 
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Assumed cross-section 




Cross-section 
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Center of mass 
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Figure 7.6: Schematic diagram of a hypothetical cross-section evolution with centre 
of mass energy and the interpolation used to evolve a cross-section at the highest 
energy down to lower energies. If the real cross-section is larger at the highest 
energy it is also higher at all lower energies. 

Two scans over the (m x , m^, tj) space were performed, with neutralino masses 
ranging from 68 to 103 GeV, slepton masses ranging from 67 to (m x — 1) GeV and 
values of log 10 (7]-) (where t\ is measured in seconds) ranging from —11.25 to —6.75. 
The step size was 1 GeV in both the masses (i.e. integer values were used). The 
values of log 10 (r^-) were non-equidistant. For one scan the neutralino branching 
ratios to each slepton were set equal at ~ (i.e. slepton co-NLSP scenario), for 
the other the neutralino was assumed to decay exclusively to stau plus tau (i.e. 
stau-NLSP scenario). The energy binning for the data was as shown in Table |7.1| 
(i.e. seven bins from 188.6 to 206.7 GeV). At each point the partial efficiencies 
were calculated for each channel, energy and optimisation using the interpolation 



and extrapolation techniques described in Section \7.4\ . The total efficiencies were 
reduced by one standard deviation (and the partial efficiencies scaled accordingly) 
to be conservative. The background expectation values (the b Se of Eq. |T7^) were 
linearly interpolated from the values determined from background Monte Carlo 
data at 189 and 208 GeV. Note that in the case of an event being observed at an 
energy and with a selection set for which 6 Se = 0, the value of the test statistic 
is infinite. This is clearly not sensible since the finite number of Monte Carlo 
background events analysed means that although no events may have been observed 
with a given selection set at a given energy, it does not necessarily mean that 
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the mean expected number is non-zero. Thus, since approximately 10 times the 
number of events expected in data were analysed at 189 GeV, and 100 times at 
208 GeV, a minimum value was placed on the bs e of (1 — t)/10 + t/100, where 
t = (\/s e — 189)/(208 — 189) and y/s e is the energy of the respective bin. A 
trial value for the cross-section limit at {y/s max =) 208 GeV was run down to the 
individual values at the seven bin energies. This gave all the information required to 
calculate a test statistic, as defined by Eq. |7^, for any given experiment outcome. 

Since there are two optimisations, there are two confidence levels. The confi- 
dence level chosen is the one which gives the higher average expected value under 
the background-only hypothesis (with no reference to the actual search outcomes). 
The exact procedure is as follows. Firstly the two test statistics for the actual 
observed outcomes under both optimisations are calculated. Then a series of 100 
'toy' Monte Carlo experiment outcomes are generated under the background-only 
hypothesis for each optimisation (i.e. using random numbers, 100 sets of hypothet- 
ical ns e are generated in line with the values given by the Monte Carlo background 
expectation), and the test statistic is calculated for each. Then a sequence of toy 
outcomes are generated in the same way but under the signal-only hypothesis, and 
the test statistic is calculated for each. For each successive pair of signal outcomes 
the confidence levels for the two actual observed outcomes and the two sets of 100 
toy background outcomes are updated. The confidence level for a given actual or 
toy background outcome after N signal-only outcomes have been generated is given 
by 

_. Tllower 
CL = 1 " IT ' 

where ni ower is the number of the toy signal outcomes which gave a test statistic 
less than or equal to the test statistic of the outcome under consideration. This 
obviously has an associated binomial error. Once the average of the 100 confidence 
levels under one optimisation becomes larger than that of the other's to a good de- 
gree of statistical significance, the optimisation corresponding to the lower average 
is dropped and the process continues with just one optimisation. The upper and 
lower limits for the value of the confidence level corresponding to the one actual ob- 
served outcome still under consideration are taken as the upper and lower bounds 
of its 99.9% central confidence interval: CL high and CLi ow respectively. These are 
calculated for N = 30, 100, 300, 1,000, 3,000, and 10,000. If CL high is less than 
95% then the point is considered un-excluded and the procedure stops. If CLi ow is 
greater than 95% then the point is considered excluded and the procedure stops. If 
neither of these is true the procedure continues until the next value of N at which 
they will be re-calculated. If N reaches 10,000 with neither of these conditions 
being met then it is not known whether the confidence level is above or below 95%, 
only that it must be close. This is the requirement that must be satisfied for the 
trial value of the cross-section at 208 GeV to be considered the cross-section limit. 
Figures [T7| to [TP] show the cross-section limit (cr 95 ) at \fs = 208 GeV in units of 



nanobarns as a function of neutralino and slepton mass for each value of log 10 (r[) 
considered. The left and right plots correspond to the slepton co-NLSP and stau- 
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NLSP scenarios respectively. The contours show lines of constant cross-section 
limit. The x (m x ) and y (m ; ~) scales are identical (and so the line of slepton- 
neutralino degeneracy is the 45° line joining the corners), and the cross-section 
limit is calculated at the mass-bin centres. The colour shading scheme is the 
same for all the plots, and the numbers on the keys to the right of each refer to 
values of log 10 (cr 9 5) (where <7 95 is in nb) and range from 0.924 (a 95 = 8.4 nb) to 4.1 
((T95 = 12, 590 nb). In the lowest and highest lifetime plots for each scenario no limit 
exists beyond the outer contour (at 095 = 10, 000 nb). The trial cross-section values 
were required to converge on 095 within a maximum number of iterations in order 
to ensure CPU-time was kept reasonable. In a few bins this did not occur, and so 
there are some spurious points where no limit is shown (appearing as isolated black 
squares). These are purely random occurrences and not the result of any physical 
effect. 

The region of maximum sensitivity moves from bottom left to top right with 
increasing lifetime, since when the lifetime is short (long), low (high) neutralino 
and slepton masses are preferred to increase (reduce) the slepton decay length 
such that it is in the most sensitive region. Figure [7.10| shows the confidence level 
distributions for each scan. They are centred on 0.95 (95%) with a small spread, as 
expected. For comparison with the quoted limits, a perfect search (100% efficiency 
and zero observed events) would yield a limit of 



a 95 = -hi(l-0.95) 



(208 - 2m x ) 
Yu^eiVSe - 2m 



X) 



where the 0.95 is the confidence level, and elements in the sum over energy bins 
for which y/s e < 2m x are ignored. The m x dependence comes from the assumed 
linear evolution of the cross-section with y/s. This function is plotted in Figure [7.11 



together with a comparison with the obtained limit at logTj- = —8.75 under the 
co-NLSP scenario. 
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Figure 7.7: Model independent limits on the neutralino production cross-section 
at a LEP energy of 208 GeV as a function of neutralino and slepton mass. The 
limits are at 95% confidence level and in units of nanobarns. 'Co-NLSP' means the 
neutralino was assumed to decay to the three sleptons with equal branching ratios 
of 3. 'Stau-NLSP' means the neutralino was assumed to decay exclusively to the 
stau. The slepton lifetime is the V in log 10 r and is measured in seconds. 
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Figure 7.8: Model independent limits on the neutralino production cross-section at 
a LEP energy of 208 GeV as a function of neutralino and slepton mass. Limits are 
given at 95% confidence level in units of nanobarns for the co-NLSP and stau-NLSP 
scenarios and a range of slepton lifetimes. 
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Figure 7.9: Model independent limits on the neutralino production cross-section at 
a LEP energy of 208 GeV as a function of neutralino and slepton mass. Limits are 
given at 95% confidence level in units of nanobarns for the co-NLSP and stau-NLSP 
scenarios and a range of slepton lifetimes. 



7.6 Model specific parameter space exclusion 



A scan over the parameter space of the model described in |2(J was performed 
by Aran Garcia-Bellido in order to translate the results of all the ALEPH GMSB 
searches into excluded regions of the model's parameter space. This was an update 



of the scan described in |L6[ using the new highest energy LEP data and including 
the analysis described in this thesis. The scan ranges for the parameters (which 
are described in Section |2.6.4|) and the number of points used in each are given in 



Table |7.3| . The five values for the upper limit on A correspond to increasing values 
of N5, i.e. for N5 = 1 the upper limit on A is 75TeV, for N5 = 2 the upper limit 
is 60TeV, etc. At each point the ISAJET 7.51 program |48|] was used to calculate 
the sparticle masses, branching ratios and lifetimes, together with their production 
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Confidence Level Confidence Leve 

Figure 7.10: The distributions of the confidence levels by which the cross-section 
limits are excluded. They are centred on 0.95 with a small spread. 
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a). b). 

Figure 7.11: a). The cross-section limit at 208 GeV, as a function of m x , that 
would result from a search with 100% efficiency and which gave zero observed 
events. The limit increases smoothly with m x because, for a given cross-section 
at 208 GeV, a higher neutralino mass causes a lower cross-section at the lower 
energies (according to the linear evolution scheme described earlier). The gradient 
has discontinuities however, corresponding to the neutralino mass exceeding the 
beam energies of successive energy bins. b). A plot of the ratio of the limit for the 
co-NLSP scenario with a log?]- = —8.75 to the ideal limit shown in the plot on the 
left. It shows that the limit obtained is close to the theoretical minimum. 
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Parameter 


Lower limit 


Upper limit 


No. of points 


M 

lvJ -mess 


1U Kjii V 


1 PI 12 PpV 


o 


Mg 


10- 9 - 5 GeV 


10- 4 GeV 


12 


A 


2TeV 


75, 60, 45, 45, 45 TeV 


151 


tan p 


1.5 


40 


9 


N5 


1 


5 


5 


sign(/i) 




+ 


2 



Table 7.3: The ranges for the scan parameters and the number of points in each. 
815,400 points were considered in total. 



cross-sections at the seven LEP energies given in Table 7.1 



These were then passed to the same routine that calculated the confidence level 
for the cross-section limit scan, and the same method was used for judging whether 
the confidence level was greater than, less than, or not distinguishable from 95%. 
Points with a confidence level greater than 95% were considered excluded, points 
with a confidence level less than, or indistinguishable from 95% were considered 
unexcluded. 

One problem arose in this scan which did not arise in the cross-section limit 
scan however. In the cross-section limit scan only the case of perfect slepton de- 
generacy, and the case of a stau so much lighter than the selectron and smuon 
it was effectively the sole NLSP were considered. In this scan, points in-between 
these two extremes exist, where the stau is significantly lighter but the neutralino 
branching ratios to the selectron and smuon are still quite significant. This does 
not pose a problem for the ee, jljl, ff and efx channels since the two sleptons in- 
volved in each channel are still degenerate with each other. But events in the ef 
and (If channel will involve sleptons of different properties in the same event - a 
case for which signal Monte Carlo data has not been generated. However, since the 
efficiencies (partial and total) depend almost solely on the slepton decay length, 
it can be expected that the true efficiencies will lie somewhere between those cal- 
culated based on both sleptons having the stau's properties, and those based on 
both having the other slepton's properties. Although it cannot be expected that 
the true efficiencies can be obtained by any kind of simple interpolation between 
the two, it can be expected that the true confidence level will lie between the 
confidence levels calculated assuming both sleptons have the one set of properties 
and both having the other. So two confidence level calculations are run in parallel 
for both optimisations, and the lowest is chosen as the real confidence level to be 
conservative. 

Significant regions of excluded parameter space in the (m x , m?) plane were only 
found for four of the twelve gravitino masses used in the scan, and these regions 
are shown in Figure |7.12| . The sudden rise in the lower edges around m x ~ 90 GeV 
corresponds to the threshold for neutralino production beginning to exceed the 
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energies at which data was taken. The positions of the lower edges of the regions for 
Mq = 10~ 8 , 10~ 7 ' 5 and 10 -7 GeV are predominantly defined by the lowest slepton 
mass for which Monte Carlo data was generated rather than any physical effect[|. 
For Mq = 10 -6 ' 5 GeV the lower edge is higher. This is because the slepton lifetime 
is proportional to (M^ 2 / 'ra~) . So for a high gravitino mass and decreasing slepton 
mass the slepton decay length can rise beyond the maximum value for which there 
is sensitivity before the slepton mass drops below the minimum defined by Monte 
Carlo data. A similar effect accounts for the fact that the right hand edge of the 
Mq = 10~ 8 GeV region is at a lower value of m x than for the other regions. Here 
the slepton lifetime is short, and a high neutralino mass exacerbates this to produce 
a decay length below the minimum value for which there is sensitivity. 



7.7 The results in a wider context 



The analysis described in this thesis was incorporated in an update of the common 
interpretation of all ALEPH GMSB work in This includes a description of a 
scan essentially identical to that described in Section \l.Q[ and contains plots similar 
to those in Figure 7.12 but showing the regions excluded by every search. The 
search described in this thesis did not exclude any significant amount of parameter 
space that the search for direct slepton pair production over the same energy range 
(||50||) did not, since although the neutralino cross-section was often higher, the 
slepton cross-section was never so low that points could not be excluded by the 
slepton pair-production search when they could be by this search. The search for 
pair-produced sleptons is also not limited by the neutralino mass and so covers 
a wider range of the parameter space. This search then serves to confirm the 
exclusion of the relevant region of parameter space under the model considered, 
but still could prove independently useful should models be hypothesised in the 
future which predict a greater disparity between the neutralino and slepton cross- 
sections. It also brings the overall ALEPH GMSB search effort close to complete. 

as: 



The current ALEPH limits are given in [49 



In the neutralino-NLSP scenario: 
99 GeV on m x for short lifetime 
60 GeV on m x for cr x < 100 m 

54 GeV on m x for long lifetime (indirect limit from MSSM chargino searches) 

In the stau-NLSP scenario: 
74.7 GeV on for any lifetime 
99 GeV on m v for no stau lifetime 



1 Note that this does not correspond exactly to the actual lowest slepton mass for which Monte 
Carlo data was generated. The technique for interpolation in y/s allows efficiencies to be calculated 
for masses down to ~ 189/208 x the lowest slepton mass. 
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Figure 7.12: Excluded regions of (m x , m=) space for four different gravitino masses. 



7. 7. THE RESULTS IN A WIDER CONTEXT 



111 



In general: 


54 GeV 
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the NLSP mass for any NLSP lifetime 
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LEP and the associated detectors have now been dismantled in preparation for 



the installation of the new LHC, and so [H9fl constitutes the final word on GMSB 



detection with the ALEPH detector. A paper combining the GMSB results of all 
the LEP experiments is being planned. No evidence for supersymmetry was found 
by LEP in any guise. The Tevatron, LHC and FLC present the next opportuni- 
ties for supersymmetry to make itself known, and their findings are awaited with 
anticipation. 
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Appendix A 
Signal Monte Carlo 



The exact points in (m x , m ; ~, dj) space at which signal Monte Carlo was generated 
for each LEP energy are given here. 

At 189 GeV the slepton decay lengths (in units of cm) were chosen for each 
channel as follows: 



ee: 0.2, 0.5, 2, 5, 30, 80, 150, 500. 

fifi: 0.2, 0.5, 2, 5, 30, 80, 300, 2000. 

ff: 0.2, 0.5, 2, 5, 30, 80, 300, 2000. 

e/2: 0.2, 0.5, 2, 5, 30, 80, 300, 500. 

ef: 0.2, 0.5, 2, 5, 30, 80, 300, 500. 

jlf: 0.2, 0.5, 2, 5, 30, 80, 300, 500. 



At 208 GeV the decay lengths 

0.2, 0.5, 2, 5, 30, 80, 300, 500, 
were used for all channels. 

The mass points at 189 GeV, in units of GeV and in the format m x -m z ~, were 
chosen for each channel as follows (the points in brackets were only generated at 
the two highest respective decay lengths): 



ee : 69-68 74-68 79-68 84-68 89-68 94-68 
74-73 79-73 84-73 89-73 94-73 
79-78 84-78 89-78 94-78 
84-83 89-83 94-83 
89-88 94-88 
94-93 



/2/2: 70-69 82-69 94-69 
82-81 88-81 94-81 
94-93 

(94-87 94-75) 



ff: 70-67 82-69 94-69 
82-79 88-81 94-81 
94-91 

(94-87 94-75) 



e/2 : 70-69 82-69 94-69 
82-81 88-81 94-81 
94-93 

(94-87 94-75) 
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ef : 70-67 82-67 94-67 
82-79 88-79 94-79 
94-91 

(94-85 94-73) 



jir: 70-67 82-67 94-67 
82-79 88-79 94-79 
94-91 

(94-85 94-73) 



At 208 GeV the mass points were the same amongst the three channels that did 
not involve a stau, and amongst the three that did. They were as follows: 

No stau: 70-69 87-69 103-69 Stau: 70-67 87-67 103-67 

87-86 95-86 103-86 87-84 95-84 103-84 

103-102 103-100 

(103-94 103-77) (103-92 103-75) 
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